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Abstract 

A survey devoted to Ao-condensate in gauge theories at high temperature is 
presented. Both the theoretical foundations of the spontaneously generated con- 
densate and known methods of its calculation are discussed. As most important 
consequence the SU(N) global symmetry breakdown is investigated in details. 
Influence of Aq on matter fields is studied in different aspects. Some new results 
concerning this subject are reported as well. 
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1 Introduction and Motivation 



In the last few years an interest to problems of high-temperature QCD has considerably 
increased. Among commonly discussed topics on the deconfining phase transition and 
the chiral symmetry restoration, new possible phenomena - generation of the gauge 
classical field (the so-called v4 -condensate) and spontaneous breaking of the global 
gauge symmetry caused by < A > - have become of great importance. A number 
of new essential results obtained in various approaches as well as interesting ideas 
and hopes connected with the non-zero vacuum value < A >^ impelled us to 
undertake an attempt to summarize here the most essential achievements and unsolved 
questions in this area of gauge theories. To be more transparent, we are going to 
formulate a general point of view both on physical nature and mathematical aspects 
of these phenomena representing them in different methods of calculations. We would 
like to discuss as well the most exciting problems and possible ways of their solving. 
Naturally, we are aware that it is a sufficiently difficult task to present a general survey 
concerning a particular question in a so rapidly developing area as gauge theories at 
finite temperature. Moreover, it cannot be excluded that by the moment of appearing 
our paper in the journal some questions discussed here will be solved. The situation 
gets complicated by the fact that as far as we know there is no (mathematically) 
strict proof at the present time that A -condensate must fall at high temperature. 
Nevertheless, we believe that calculation of < Aq > by different methods and, on the 
other hand, derivation of the most significant consequences of such a condensate on 
multi-particle systems make our attempt quite justified. Besides, the appearance of 
A -condensate and the breakdown of global gauge symmetry can undoubtedly lead to 
significant improvement of our conception both of the high temperature behaviour of 
the strongly interacting matter and of the physics of the gauge theories on the whole 
and surely have a connection to other problems currently under investigation (as, for 
instance, infrared problem, behaviour of quarks at non-zero baryonic number, etc.). 
All these questions will be considered in the paper. 

Let us begin with a comprehensive consideration of some known facts obtained 
from the studies of QCD. The Hamiltonian can be formally written as the sum of the 
chromoelectric and the chromomagnetic terms and has the following form in lattice 
version of the theory 

H = 7T E2 + 7^UUU+U + = H E + H B . (1) 

At finite temperature the behaviour of chromoelectric fields has been well studied 
both in the perturbative (in g 2 ) and especially in the non-perturbative regions. As is 
generally known, in the strong coupling approximation the main contribution to the 
partition function results from the chromoelectric part in eq.(^), because the chromo- 
magnetic term being proportional to g~ 2 can be treated perturbatively in jp. At high 
temperature because of periodic boundary conditions the gauge field configurations 
called Polyakov loops develop a non-vanishing expectation value, which breaks Z(N) g i 
symmetry of the initial QCD-action and leads to the deconfinement. As the Polyakov 
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loops transform non-trivially under Z(N) g i rotations, their non-zero expectation value 
could mean screening of Z(iV)-charges (or static quarks) at T > (where T C D is the 
critical temperature of the deconfmement phase transition). If this is the case, one 
may claim that there exists a physical quantity which characterizes the phenomenon 
of screening. This quantity is called the Debye mass and is defined in the continuum 
as the zero momentum limit of the zero-zero component of the vacuum polarization 
tensor, 

m 2 D {T) = -n 00 (£^0,A; = 0). (2) 

This definition gives a gauge invariant value only in the lowest non-trivial order of the 
weak-coupling expansion for SU (N c ) gauge theory with Nf massless fermions 

/V N t 

-U 00 =(f + -^)g 2 (T) + O(g 3 ). (3) 

Calculation of Eq.(^) even in the two-loop approximation reveals infrared divergences 
and gauge dependence |1J] leading to the conclusion that such a calculation cannot be 
trusted as high-order corrections are not calculable systematically. The starting point 
of the lattice studies is the connected correlation function 

T(R,T) =< W(0)W(R) > - < W(0) > 2 (4) 

of the two Polyakov loops 

N t 

W(x) = Spl[U (x,t) (5) 

t=i 

separated by a spatial distance R. The colour averaged quark- ant iquark potential V(R) 
defined by this correlation 

V(R) = -lnT(R,T) (6) 
has been found to be of the screening form [0 

con 

V(R)^——exp(-m D R) (7) 
K 

In such a way we obtain an other definition of the Debye mass. In the weak coupling 
approximation the definitions @) and (0) have to coincide |§. 

Results of lower accuracy have been obtained from studies of the chromomagnetic 
fields behaviour at high temperature. First of all, two essential problems being im- 
portant for understanding the finite-temperature physics of non-abelian fields on the 
whole are to be emphasized: infrared problem and area law for spatial Wilson loop. 
The infrared problem appears at attempts to study perturbatively the physics of the 
high temperature region: The static chromomagnetic sector develops infrared diver- 
gences which make the perturbative expansion invalid starting from g e order. Leading 
infrared divergences are those of the three-dimensional Yang-Mills theory [[|. It is 
widely believed that some physical characteristic should exist leading to screening of 
the chromomagnetic forces and to curing of the divergences. It is usually called the 
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chromomagnetic mass though there are several definitions of this quantity in the liter- 
ature. 

1) By analogy with the definition of the Debye mass (@) the same limit of the space 
components of the H^ko, k) could be considered: 

mJ(T) = -U u (k^0,k = 0) (8) 

However, this limit was found to be equal to zero in the one- loop approximation [H. 
Two-loop calculations performed in the axial gauge have given finite result (but dif- 
ferent in various schemes) ||. Besides, specific difficulties of the axial gauge (coming 
from its singular character) bring an ambiguity in obtained results and, what is more 
important, lead to qualitatively different behaviour of Uu in comparison with relativis- 
tic gauges. The most general two-loop calculations of 11^ in an arbitrary relativistic 
gauge were done in 

„f2ur , \ 7iq A N 2 r a 3 4a 2 13a . , k 2 

- nSfc*. = o) = sttwiT + - + — ~ 21 "V (9) 

where a is the gauge fixing parameter (a = 1 corresponds to the Feynman gauge). 
There is the infrared singularity in the last equation when k 2 — > 0. One can see as 
well that the sign of the right-hand side of (|j) can be changed by varying a. Hence 
one must conclude, expression @ cannot define a physical mass. A detailed discussion 
of such a situation may be found in J?|. Consequently, at the present time we do not 
possess any regular methods for calculation of the chromomagnetic mass (defined in 
the sketched way). That is why other parameters should also be discussed. 

2) By analogy with the definition of the Debye mass @-(|7|) the chromomagnetic 
mass can be determined from the corresponding correlation functions. In this case it 
has to be the spacelike Wilson loop in the adjoint representation. If m M 7^ 0, "heavy" 
gluon current sources have to be screened. As far as we know such a calculation has 
not been performed yet. 

3) If gluons acquire a magnetic mass, the magnetic fields of test charges will be 
screened. It means we can measure a monopole-antimonopole potential with the 
charges in the center of the gauge group. Corresponding Monte-Carlo simulations, 
although on the small lattices, were presented in ||. The "magnetic" potential really 
is of the screening form and the chromomagnetic mass determined from this potential 
is in a good agreement with the formula m M = Bg 2 (T)T, where g 2 (T) is the coupling 
constant renormalized at the momentum scale T. The computations were performed 
both for the Wilson pure SU(2) gluodynamics and for the Mack-Petkova modified 
SU (2)-model which does not contain dynamical monopoles. 

It seems to be a logical conclusion that all of these three definitions determine the 
same physical quantity though we do not know a strict proof of this statement. 

Second important (and strict) result which should be understood and obviously 
concerns the chromomagnetic fields was obtained by Borgs (see also |K|). Its 
gist is the following: The spacelike Wilson loop obeys the area law at arbitrary high 
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temperature. In principle, this result does not look very strange. The chromomagnetic 
sector in Eq.(|l|) at finite T does not generate (unlike the chromoelectric part) any 
configuration transforming non-trivially under Z(iV) 9 2-rotations. Therefore it seems 
natural that there is no screening effect of the corresponding sources (currents) being in 
the fundamental representation in the vacuum. Two points we would like to emphasize 
here are the following: 

1) area law explains the impossibility to apply the standard weak coupling expansion 
at least to all modes of the chromomagnetic sector because the area law can never be 
achieved in this expansion. 

2) as is known from the duality relations [IT], [[HJ, the area law for the Wilson loop 



implies the perimeter law for the 't Hooft disorder parameter. In its line this has to 
signify a screening potential for a monopole-antimonopole pair. 

Thus we can see that despite the fact that there is no strict proof that the chromo- 
magnetic mass should exist we have to take into account that some kind of screening 
of the chromomagnetic forces does exist. We consider that the result by Borgs and the 
duality relations surely lead to this conclusion. 

The crucial question coming from these facts concerns the nature of this screening. 
Undoubtedly, a mechanism of such a screening implies the existence of the magnetic 
mass in one of the meanings above. 

Two possible mechanisms supplying gauge fields with mass and leading to the 
screening of chromomagnetic forces are available and discussed now in the literature: 
dielectric structure of the QCD vacuum and the spontaneous breaking of global gauge 
symmetry in the background < Aq >. 

1) Dielectric vacuum. 

The attempts to formulate a confinement model in QCD framework led in the 
second half of the seventies to the so-called dielectric theory. In this model the gauge 
action is accompanied by additional terms corresponding to the auxiliary field p. If 
we calculate Gauss' law from the action of this kind we shall formally find the relation 
with dielectric field as in the electrodynamics of the dielectrics. The definition of the 
pure gauge model is the following |T5| 



X,fl 



Z= [l[D(<S>)eM-Sm (10) 

where 

s@) = E E + ^ E W (*) + E vw>\ < ii i 
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= P^, Up E SU(N), < p < oo (12) 

£>($) = p 3 dpDp(U), (13) 

Dp{U) is the invariant group measure. The dielectric field is introduced via the auxil- 
iary field p. V(p) is a gauge invariant potential for this field, p is a scalar field in the 
colour space for SU(2) and a tensor field for SU(N > 2). One of the most essential 
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features of such theories is a possibility to introduce a gauge-invariant mass for gauge 
fields $fj,(x) (see the second term in (|TT|)). 

It has been recently shown that the lattice effective action for the infrared dangerous 
static modes appears at the dimensional reduction to be a little complicated version of 
the described dielectric theory ||14|| . Physical mass of the dielectric field is 



m d = T(A° - 2A e cos\a p g <A >)) (14) 

where A = 2 (fiat function of T, g), A e = p^-. Undoubtedly, this mass leads to the 
screening of chromomagnetic forces and therefore of all gluon sources being in the 



adjoint representation |14j . It means in a sense the solution of the infrared problem, 
truly, perturbative expansion is invalid in any case since the dielectric field appears 
only at the non-perturbative level and vanishes if we try to apply the perturbative 
expansion. 

2) v4 -condensate. 

The second possibility is the spontaneous breaking of the SU(N) global symmetry 
at high temperature accompanied by the generation of the chromomagnetic mass. Just 
this phenomenon will be in the focus of our interest here. The only known mechanism 
providing such a breakdown at high temperature is the generation of Ao-condensate 



which is a constant part of the temporal gauge field component Aq = const ||15||. The 



nature of this parameter is connected with special properties of non-abelian gauge 
fields at finite temperature which display themselves in different ways depending on 
the particular method of calculations. We adduce below a comprehensive description 
of these properties. In general, the possibility of < A >^ originates from the 
compactification of the imaginary time direction at T ^ 0. So, it is necessary to 
combine periodicity in the imaginary time A^(0) = ^4^(0) = A' J/3) and the 

gauge transformations of the fields A'^ = UA^U + + ^Ud^U + where U is the gauge 
transformation operator. To make this transformation agree with the periodicity, the 
operator U = exp(igA /3) (where A = A^tf 1 ) must commute with generators t a . Hence, 
A should belong to the center of the gauge group: A = J^, n e Z N ,n = 0,1, N — 
1. Actually, the vacuum value of A has to be calculated from a full effective action with 
quantum fluctuations included and if it will be found that < A > differs from J^p, 
the spontaneous breaking of the gauge symmetry will be determined. In what follows, 
speaking about Ao-condensate we will assume it as a non-zero vacuum value of zero 
gauge field component obtained at a minimum position of the full effective action. The 
presence of the classical gauge field in vacuum may give a necessary missing parameter 
for solving the infrared problem. Moreover, this condensate effects various processes 
at high temperatures. In particular, it may lead to the spontaneous breaking of the 
baryon charge symmetry, etc. 

As a matter of fact, the two screening mechanisms could be somehow connected. So 



it has been discussed in that the mass of the dielectric field in the reduced theory 
can be proportional to < Aq >. We shall consider this possibility in the corresponding 
part of the survey. 
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The described general idea can get different realization depending on applied cal- 
culation methods. A general picture of the effects of the A -condensate will appear 
when the results obtained in various particular approaches will be gathered together. 
The goal of the present review is to describe in a systematic way the results obtained 
in three approaches to A -condensate calculations. We discuss the method of effective 
Lagrangians for A -field, the Hamiltonian description of the SU (iV)-symmetry break- 
ing in the background A both on the lattice and in the continuum theory, and the 
loop expansion of the effective action in QCD. In all these approaches the gauge field 
condensate has been determined in the framework of the appropriate approximation 
schemes. Anyway, at the present time there is no common opinion about realization 
of this phenomenon in the nature because of a number of problems concerning the 
accuracy, gauge invariance and precision of the calculations. So, in what follows we are 
going to adduce the comprehensive analysis of these problems and compare the results 
obtained. 

Our Survey is organized as follows. 

First we remind in brief the most essential features of gauge theories at finite tem- 
perature and introduce our notations (chapter 2). The general status of A -condensate 
is presented in chapter 3. In chapter 4 we go through the calculation of the conden- 
sate in the Hamiltonian formulation on the lattice. The gauge independence of the 
condensate in this approach is discussed as well. In chapter 5 the loop expansion for 
< A > calculation is examined. The central point of this examination is the proof of 
gauge independence of the condensate by means of the Nielsen identities. Chapter 6 
is devoted to elaborating of the effective Lagrangian method for the expectation value 
of Aq. In chapter 7 we compute the condensate in the theory with dynamical quarks. 
We consider both lattice and loop approaches and compare their results. Some conse- 
quences of the non-zero < A > are reviewed in chapter 8. Here we are going to give 
some new results related to the A -condensate phenomenon. So we give a sketch of 
the calculations of the heavy quark potential in the background < Aq >. The infrared 
problem is reexamined at < Aq >^ in reduced lattice theory at high temperature. 
We give a proof that the adjoint spatial Wilson loop obeys perimeter law, which means 
the screening of chromomagnetic forces. Further, we discuss the spontaneous breaking 
of the baryon charge symmetry and some closely related problems. The general con- 
ditions for appearing of the Chern-Simons action in the background < Aq > at high 
temperature are presented as well. Brief Summary and Discussion can be found in 
chapter 9. 

2 Gauge theories at finite temperature. General 
outlook 

In order to make our survey as self-contained as possible and for convenience of readers 
we would like to make a sketch of SU (N) gauge theories at T ^ 0. To begin with, we 
remind the general formulation of the finite temperature theories both in the continuum 
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and on the lattice and introduce our notations. In the continuum the QCD action has 
the form: 



S = J[d 3 x]J^ dtLgcD, 



Lqcd = —TrF^F^ + J2^ f [Dnln-m f + (D + lt i) l0 ]^^ (15) 
z 9 /= i 



where 



= dfj,A u - d v Ap + ig[A^ A v ], D^ = d^~ igA„ A„ = Ap\ (16) 

H is the baryonic chemical potential. 

On the lattice we consider the Wilson action: 

s = s G + s f = 

2 ^~JrU { a P)] + aa + 

x,n=-d f x 

where U n (x) G SU(N) and Tr U(dP) is the plaquette character in the fundamental 
representation. The form of the matrices T XjX+n and T depends on a sort of lattice 
fermions. 

It is well-known that the problem of quantization of gauge theories depends essen- 
tially on the boundary conditions imposed on the gauge fields. The compact topology 
leads to another physical picture of the gauge model than the Euclidean topology 
accompanied by the corresponding boundary conditions. QCD at zero temperature 
is studied in the Euclidean space G = R d . A transition to the cylinder topology is 
achieved by imposing the periodic boundary conditions on the gauge fields along the 
imaginary time direction with period (3 = l/T which plays the role of inverse temper- 
ature. 

Thus, the finite temperature theories are defined on the space G = R^ 1 <g> S 1 . The 
periodic boundary conditions generate the new observable known as the Polyakov loop 
0: 

rP 

W = P exp (ig / dt Ao{x,t)) (18) 
Jo 

(or see @ for the lattice definition). 

The finite temperature formalism is used for studying the thermodynamical features 
of QCD and its different phases. So, let us discuss in brief these features and, first of 
all, the phase structure obtained from lattice calculations. 

Nowadays two phase transitions are usually considered to be relevant to QCD at 
finite temperature and/or baryon density. One is the deconfinement phase transition, 
i.e. transition from the confining phase of the hadronic matter to the quark-gluon 
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plasma phase. The other is the chiral symmetry restoration phase transition. As is 
generally believed, the physical picture obtained from lattice studies could be the fol- 
lowing. The pure gluonic QCD action has the global Z(N) symmetry, which leads to 
the confinement of color quark charges in the low temperature phase. In the high tem- 
perature phase Z(N) symmetry is spontaneously broken and quarks are screened by 
some specific configurations of the gluon fields. This transition is called the deconfine- 
ment phase transition and it is best studied in pure gluodynamics. Yaffe and Svetitsky 
T7| have proved the relation between the free energy of the infinitely heavy quark F q 



and the expectation value of the trace of the Polyakov loop, < TrW >= exp(— f3F q ). 
This relation and the behaviour of TrW x under Z(N) g i transformations 

TrW x — ► zTrW x (19) 

were used as the basic ones to study the deconfinement phase transition in the Yang- 
Mills theory [TS]. Difference of the expectation value of the TrW x from zero is the 



signal of the spontaneous breaking of Z(N) g i symmetry and also of the deconfining of 
the static colour charges transforming non-trivially under Z(N): 

t i n-^TtW \\ - J °' T < T ?> confinement phase, , . 

{i xi) 1 Z * f(T), T > T C D , Z e Z(N) gh f(T) < 1, deconf. phase 

In the full QCD with dynamical quarks Z(N) global symmetry is evidently broken 
from the beginning and no other appropriate order parameters for this transition are 
known. On the other hand, in the presence of the dynamical fermions the first order 
deconfining transition in the pure SU (3) gauge theory is getting weaker as the quark 
mass decreases. This weakening is going on up to the quark mass of order T. However, 
for smaller masses the transition has been found to become stronger with the quark 
mass decreasing. This transition can be characterized by the chiral order parameter 
< a > and is called the chiral phase transition. Numerical results show that the decon- 
finement and the restoration of the chiral symmetry, being distinctive by character at 
the first sight, are in fact indistinguishable in the region of intermediate quark masses. 

Before proceeding further to discuss the gluon field condensation at high temper- 
ature let us make a short summary on the perturbative QCD vacuum at T ^ (see 
fT5|| , fL9"| , [p0|1 for a detailed review of the finite temperature properties of QCD) The 
perturbative vacuum is supposed to have no structures or condensates. Its main prop- 
erty is the asymptotic freedom (g{T) —>■ in the limit T —>■ oo |15[). At T ^ the 
momentum space is naturally divided in two parts: 1) | k |^> gT where the perturba- 
tive methods and results are to be reliable owing to asymptotic freedom; 2)| k |<g[ gT , 
this is a truly infrared region. Just for these momenta one runs into infrared divergen- 
cies in the higher orders of perturbation theory, gauge dependence and other problems 
which signal that the perturbative vacuum is not adequate to the nature. So, as is 
expected and has been discussed in the introduction, some new macroscopic parameter 
(like the Debye mass) should be dynamically generated. As a candidate, the gluon 
magnetic mass, calculated by perturbative and non-perturbative methods [2I|, |J has 
been discussed. 
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At the same time another possibility has also originated from the analysis by per- 
turbative methods. Here we have in mind the dimension reduction process at high 
15fl , [f2~2]l , p3f , ||24j| . The basic reason is as follows: in the gluon free 



temperature 



propagator k 2 + ^TmT) 2 °f imaginary-time perturbation theory, the term limT acts like 
a "mass" in the three-dimensional theory. According to the Appelquist-Carazzone the- 
orem p2| all nonstatic modes n^O decouple in the limit T — + oo, leaving the static (= 
three-dimensional) sector as the effective theory. This idea can be realized in different 
ways and, actually, a number of problems of calculation of the effective Lagrangian 
appeared [24 1 . But in general it gives a possibility to search for a non-trivial vacuum 
at finite temperature. In the present paper the problem of the investigation of the 
gauge field vacuum via the described mechanism will be one of the discussed topics 
both in continuum QCD and beyond perturbative horizon on the lattice. 



3 General status of Ag-condensate 

In this chapter we would like to give general reasons for condensation of the electro- 
static potential and derive the corresponding mathematical foundation. We formulate 
a mathematical task of calculation of < Aq > and give an outlook on the previous 
investigations. 

Firstly, we consider the lattice gauge theory without periodic boundary conditions 
(QCD at T = 0) with the partition function of the form: 

Z = J D/2(U)D^D^exp(-S G - S F ). (21) 

The gluonic S G and the fermionic S F parts of the action are expressed in eq.(|I7]). The 
gauge field integration is performed over SU (N)-mvariaiit group measure Dfi(U). Due 
to this fact the results of the calculations are gauge independent. For example, 

(B) = (B(U)) Uo=1 (a =o), 



if B{U) is a gauge invariant function||12||. This means that there is a gauge transfor- 
mation which allows us to fix gauge Uq = 1 or Aq = 0. We can easily see that any 
constant (Aq) may be singled out by an appropriate gauge transformation. 

The partition function at finite temperatures is calculated within the following 
periodic boundary conditions: 

U^(x,t) = U^(x,t + Np), 

*(x,r) = -VfaT + Np), (22) 
*(x,r) = -Vfar + Np) 

where Np is a number of the lattice sites in the time direction. These boundary 
conditions are incompatible with the gauge Aq = 0. We may fix a static diagonal gauge 
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only. Let Uq(x, t) = V x , where V x is a diagonal time-independent SU(N) matrix. After 
performing the gauge transformations 

U n (x,r) — (V x ) T U n (x,r)(V x+n )- T 
U+(x,t) — > (K + nrf/„ + (x,r)(K)^ 

?2 = 1, d 

*(x,r) — > (\4)^(x,r) 
%(x,t) -^V(x,t)(V x )- t (23) 

(which lead to the gauge Aq = at zero temperature) we can remove the matrices V(x) 
in the action from all links but the last one where all of them are grouped gathering in 
the Polyakov loops W x . On the lattice in the static gauge the Polyakov loop becomes 

N 

W X = P]J U (x, t) = exp(i[3gAo(x)). (24) 

T=X 

In such a way we obtain the new terms in the chromoelectric and in the fermionic parts 
of the action 

9 x,n 

S F (W) = J2%(Np-l)r x>x ± -W x (±0)* x (l), (25) 

x;±0 

W,(±0) = { ^ 

These terms are absent in the theory without periodic boundary conditions and describe 
the interactions between gauge fields, fermionic fields and the Polyakov loops. Owing 
just to these new terms the SU(N) gauge theory has the non-trivial phase structure 
described above. Moreover, we can deduce from the last equations a possibility of the 
spontaneous breaking of the SU(N) global gauge symmetry caused by condensation 
of the chromoelectric potential. The symmetry group of the pure gauge theory is the 
direct product of the group of the local gauge symmetry, the group of the global gauge 
symmetry and its center subgroup: SU(N)i oc x SU(N) g i x Z(N) g i. The last one acts 
only on the Polyakov loops. The fermionic part of the action ( p5|) violates the global 
center symmetry explicitly, the symmetry group becomes SU(N)i oc x SU(N) g i and 



the expectation value of the TrW x is not the appropriate order parameter at all|25 
Nevertheless the eigenvalues of the Polyakov loop 

W x = eMm(x)) (26) 

can be used for studying the behaviour of the system under SU (3) g i symmetry because 
the Polyakov loop transforms under 577(3) rotations like the matter fields in the adjoint 
represent at ion 

W x — >U W x U + , U e SU(N) (27) 
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The eigenvalues of the Polyakov loop can be chosen, for example, in the SU(3) theory 



as 

(/%rVi = ^o+^o, 

{Pg)- x <P* = - A ° + ^ A *> 

N 

Ew = o 

1=3 

Therefore the expectation value of W l x may be different from zero only if the spon- 
taneous breaking of the SU(3) g i symmetry is possible. The constant field A can be 
cancelled in the action only when W((Aq)) belongs to Z(N), as one may easily see 
from eq. (fffi) because in this case the matrices W commute with U n . In any other 
case the constant Aq will be present in the action and the global 5'f/(3)-symmetry is 
broken up to its Cartan subgroup 

W{(A )) 3 Z{N) =*► SU{N) — > [^(l)]^ 1 (28) 

It is clear now that such a mechanism can work both in the pure gluodynamics and in 
the theory with dynamical quarks unlike the mechanism of the spontaneous breaking 
of the Z(N) g i symmetry. 

Thus, we have two possibilities to quantize the gauge theory in the space G = 

1) The theory in G = R d ~ x ® S 1 must always possess the same symmetry of the 
vacuum and the same Lagrangian as the theory in the Euclidean space G = R d . So, for 
instance, any constant Aq should be singled out by the corresponding gauge transfor- 



mation from eq. (|23|) . This constraint leads to the restriction W l x = exp(^qi(x)),q 



0, 1, N — 1 and for the quantum theory we obtain 

r _ 2,711 - 

Z=J2 / Dfi(U)D^D^e^(-S(e^(— qi (x)),U,^)). (29) 

q t (x) JU n (T) = U n (T+P) N 

2) The global symmetries of the vacuum are determined by dynamics of the gauge 
system itself. Then, in eq.(|^) we have J2 qi ( x ) ~^ I dfi(W x ), where dfi is the invariant 
SU(N) measure. 

In the former case there is the only possibility of the spontaneous breaking of the 
global Z(N) symmetry whereas in the latter one the spontaneous breaking of both 
Z(N) and SU(N) global symmetry can occur. It is difficult to give a preference to 
any of these formulations from the theoretical point of view. The second formulation 
commonly seems to be more relevant. In what follows just this formulation will be 
explored here. 

Thus, from the picture formulated above we come to a task of calculating {A ) which 
can be formulated in the following way: The configurations A n — (n — 1, d), Aq = 
const are the solutions of the Yang-Mills equations with the additional condition that 
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the classical action equals zero. At the classical level, however, the breakdown of the 
symmetry is absent and the system is in the symmetry invariant minimum. Can the 
quantum fluctuations generating the effective potential for A lead to spontaneous 
breaking of the global symmetry and produce the minima of the effective action with 
a non-trivial value of (PgA ) ^ ^ ? On the lattice the situation is the same: for the 
configurations U n = 1, W\ = const, the classical lattice action S G is equal to zero and 
therefore we can ask the same question. 

Actually, the investigation of the A -condensate started ten years ago and < A >^ 
was determined in various approaches. Nevertheless, at the moment we have neither a 
strict proof that the condensate has to appear nor a common opinion about its genera- 
tion. This is mainly due to the mathematical difficulties which have been encountered in 
the used approaches. For example, in the loop expansion method < Aq >^ is derived 
from the two loop effective action W(A , f ) @ ,g7J, J2B|, |29| , g|, | 



1,0,11,111. 

So, to verify the result the three-loop contribution should be evaluated. This very 
complicated task has not been solved yet. Other approaches discussed in the literature 
HI,B HI, BUI also contain either some uncertainties or unsolved problems. 

Recently, the possibility of Aq condensation has been called in question from the 
point of view of its gauge invariancep9||. As was found in the background gauge, 



both the effective action and its minimum value appear to be dependent on the gauge 



fixing parameter £ 



Hence, a doubt about gauge invariance of the phe- 
nomenon has arisen. To resolve this doubt in the perturbation theory several methods 
for gauge invariant calculations have been conjectured |5fJ,|5!|,[[flJ. But in all of them 
the result - no real condensation at two loop level - has been stated. On the other hand, 
the gauge invariant results of the analytical lattice investigations |36[ unambiguously 
show that condensate does appear. The same conclusion has been supported via the 
Nielsen identities method [42| at the two loop level (in order g 2 in coupling constant). 
Another approach to the problem of gauge invariance of the condensate was recently 
proposed in]|0| where the authors built a partition function for the eigenvalues of the 
Polyakov loop in the Hamiltonian formulation of the continuum theory. The result 
A = was obtained. We shall discuss this result later on remarking that this con- 
clusion is in obvious contradiction with apparently gauge invariant lattice calculations 

The situation needs to be clarified in any way because Ao-condensate, if it does 
realize in the nature at high temperature, would be a very essential element of the 
self-consistent finite temperature gauge theory. It seems to us that the only way to 
comprehend the situation is to fix the strong results of various approaches and to 
find out both their common point and sources of discrepancies. At the same time it 
would be very desirable to find crucial phenomena connected to < Aq >^ which 
can be detected in future experiments on the heavy-ion collisions. These are main two 
purposes of the present paper. We believe that being gathered together, the results of 
various approaches can help to elucidate a lot of difficulties and outline the prospects 
for future investigations. 



13 



4 Hamiltonian formulation of the gauge theories 
and Ag-condensate on the lattice 

The first approach we would like to analyze here is the Hamiltonian formulation of the 
gauge theories. It is very desirable to have an apparently gauge invariant approach 
to defining and calculating < Aq >-condensate. Possessing this essential property, 
the Hamiltonian approach allows us to be convinced that A Q condensation will be a 
gauge invariant phenomenon. On the other hand, we want to have a strict proof that 
< A in the framework of the reliable approximation scheme. Fortunately, such 
an approximation does exist. It is the strong coupling expansion in the lattice theories. 
In what follows we consider, at first, the strong coupling region of the Hamiltonian 
lattice formulation. Then, the continuum version of the Hamiltonian approach will be 
discussed. 

The detailed description of the lattice Hamiltonian formulation can be found in the 
series of the papers |[43|| , MM and in the review fH)|| . The Hamiltonian approach on the 



lattice was developed for the first time in [46]. In the construction of the Hamiltonian 



partition function we follow our own method [|^J adducing a proof of its equivalent to 
earlier methods. The Hamiltonian of the lattice gluodynamics in the strong coupling 
approximation includes only the chromoelectric part 

H = E (30) 

links 

where E(l) = id/d(Ai) - are the chromoelectric field operators. In this approach the 
chromomagnetic term can be treated perturbatively at g 2 — > oo. Calculation of the 
partition function 

Z = Spexp(-f3H) (31) 

is connected with summing over local gauge-invariant states. This is reflected in sign 
Sp in (|31|). The corresponding physical Hilbert space is determined by Gauss' law. To 
satisfy the latter the usual method is to introduce the necessary 5-function. In order to 
solve this task we use the more general method connected with the projection operator 
technique. It should be emphasized at once that the Hamiltonian formulation allows to 
obtain an effective action for the eigenvalues of the Polyakov loop in the gauge A Q = 0. 
In what follows we are going to demonstrate that carrying out the projection onto local 
gauge invariant states and summing over the eigenvalues of the diagonal operators of 
the gauge group are equivalent to integration over spatial gauge fields in the Euclidean 
version of the lattice theory without fixing any gauge. It is a rather important point 
which was missed, for instance, in Ref. [4C]. 

Our starting point is the partition function ([JT]) where we implant a projection 
operator in each lattice site 

Z = 5pex P (-/3 £ (f )E 2 (l))P,P = Y[P° X (32) 

links Aa x 
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where the operator P r is defined as 



P r = Xt^r(9)Uj(9). (33) 
I 9 I 3=1 

To make the next mathematical construction more transparent we have put down here 
P r for a discrete group, h r is the dimension of the representation r, | g | is the rang 
of the group and Q r is the character of the irreducible representation r. The matrices 
Uj(g) are located on the links. In eq. (|3~3"D we must consider all representations which 
can be combined in such a way to form singlets in each lattice site. Eq. ( p3|) can be 
rewritten to the form 

h n 9k h n 

pr = rhH «t ,r E u kd = ArE «r • (34) 
i y i fe=i 3=1 i y i fc=i 

The second identity is the definition of the class operator C(k) which contains 
elements, n is the number of linearly independent class operators of the representation 
group G [Q. There exists a representation for C{k) 



C(k) = J %f:U J U k U. 1 (35) 

I I 3=1 

from which it is easy to deduce the most important property of the class operator: it 
is invariant under transformations out of group of representation G or in other words 

[C(k),U(g)} = 0. (36) 

A generalization of ([34]) and (|35|) on the Lie group is not very complex and is founded 
on using the so-called "unitary Weyl's trick". As is known, any unitary representation 
of the Lie group may be presented in the form 



U^) = V(v)a<P)V-\v) (37) 

by the appropriate choice of the unitary coordinate system. In these terms the invariant 
measure is 

rf$ = dvd(i(<p), (38) 
which allows to integrate over //-variables. After that we have for P r 

p r = h r J dW' r ($)£/($) = h r J dn((t))n*> r J dvU{<j>, v). (39) 

For class operator we obtain in the case of the Lie group 

C(<f>) = f duU((/),u). (40) 



It is obvious that C(0) possesses the property of the invariance analogous with 
Just this property is the most important one in the selection of the local gauge-invariant 
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states. Let C/ be the eigenvalues of the Hamiltonian ( p0|) and U^(x) be the eigenfunc- 
tions of E 2 . They are taken as variables on the lattice links, x is the lattice site and \x 
is the unit vector. Then, the partition function (|3^) is rewritten in the form 



Z = J2K l exp(- 1 C l ), (41) 
l 

where 7 = I 2 - . The operator selects gauge-invariant states among all eigenstates 
of the Hamiltonian and has the form 

Ki = J Mu)UKi p Ii u ^- ( 42 ) 



x.v 



From above formulae it follows that 

p uKu = I nwn^w^w- ( 43 ) 



x,u 



The operator P satisfies the relation PP = P and is nothing but the projection opera- 
tor. How it is working may be seen from the last three equations. Next, we would like 
to describe how to connect this technique with standard projection and with Gauss' 
law. 

Let us choose the transformations Q(4> x ) in the form 

n{<f> x )U Xt „ = exp(iE^4> a x )U XjU , 
U x ^{<j> x+V ) = U x>u exp(-iE% +1/> J a x ), (44) 

where the generators E x v are the chromoelectric field tension operators defined on the 
lattice links. Then, grouping E xv at the same (f> x we obtain that 

P = II p l P * = I «ap("M (45) 

x J 

where 

Q a x = J2 K,u (46) 

V 

are the colour charge operators and the generators of the gauge transformations in the 
lattice sites. P° in Eq. fl45|) is the projection operator onto states with zero charge and 
is a particular choice of the operator (|39|) . 

The connection of Gauss' law 5(AE) with a projection operator technique is not a 
trivial task as the calculations are quite complicated because the operators Q a do not 
commute with each other and, consequently, it is impossible to find eigenfunctions for 
all Q a at the same time. Applying "Weyl's trick" again we present S(Q a ) in the form 

. N-l 

S(Q a ) = / d^-^vd^-'^V^expli <b a Q a ]V + (v). 
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In the last equation the matrices V{y) diagonalize exp [i Ea' 1 a Q1 in the sense that 
in the final expression only Q a belonging to Cartan subgroup enter. After integration 
over v one finds 

6(Q a ) = [ dfi(<f>) exp[* Y rQ a ] = P°- (47) 

J a 

The last step we need to do is to introduce a summing over eigenvalues of the Q a of 
the Cartan subgroup. Just eigenf unctions of these operators and eigenf unctions of the 
Hamiltonian produce a full set of the local gauge invariant states in such a treatment. 

The operator P r at r 7^ projects onto the invariant states in the presence of the 
background probe charges. To prove it we have to consider the 5-function 5(Q a — q a ) 
where q a are the probe charges and the generators of the r-th irreducible representation. 
Summing over all charges belonging to r-representation we easily find 

. N-l 

Y <K<9° - <f ) =P T = / a>(0)fi*' r (0) exp[i Y <P a Q a \- (48) 

q a £r a 

Now let us consider the Wilson lattice action in the same strong coupling approxi- 
mation (restricting ourselves to chromoelectric part of the full action in flTT| ) at finite 
temperature). It is known that after integration over space gauge fields and taking a 
limit ap —>■ in the time direction we come to the same expression for the partition 
function as it appears in the Hamiltonian formulation. Since in the last case we do 
not presuppose any gauge and start from an apparently gauge invariant formulation 
let us demonstrate it once more. To integrate out the space gauge field configurations 
we follow the standard procedure and expand the Wilson action into series over the 
characters of the irreducible representations of SU (iV)-gauge group. Then, performing 
the invariant integration over dU n (x) we come to the following expression 

z G = / nw^tfflE^^^n^^'ii^nw+n.*))] ^ 

J x,t x,n r t=l t=l 

where (3 a = We can now choose Weyl's representation Q57D for the Polyakov 

loop W r = ]7t=i UQ(x,t)) rewriting the invariant measure in the corresponding form. 
After cancellation of the non-diagonal parts of the Polyakov loops we have 

Z G = [ IJ^Kt) WL C r"W{SpW r {a x )){SpW^ r {a x+n ))\ (50) 

x,t x,n r 

where W x = exp[i Yn=i a x(t)T a ], T a are diagonal generators of SU(N). It is enough 
to integrate in the last equation only over one chosen a(t'). Because of the invariant 
integration the result will be independent of a(t 7^ £'). Thus, we obtain 

z G = J n^^ntE^^^w^K+n)]- (si) 
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We consider now the quark contribution to the partition function. To get it in 
the Hamiltonian formulation we have to pick up the operator P r (|39p when r / 
and take into account both the states without charges and the contribution of the 
antiparticles when summing over q a is performed in ( |4"§|) . In such a way one finds the 
quark contribution to be 

Y[(l + Ren r x ). (52) 

X 

The same contribution can be obtained in the Euclidean formulation where we do not 
fix any gauge. In this case we consider the interaction of matter fields with Uq(x, t) and 
explore the Kogut-Susskind massless fermions Calculating the partition function 



Z q = [Y[(M x( M x ex P (-q x D q y ), (53) 

J X 

where 

A) = -{U (x,t)5y )X+0 - U(l(x,t)Sy tX -o) 
and Uq belongs to the fundamental representation of SU(N), we get 

Z q = Y[(l + ReQ x ). (54) 

X 

This contribution is exactly the same as in the Hamiltonian formulation (^) if we 
choose r e fundamental representation of SU(N). To do the last step in our proof we 
notice that at ap — > 0, Np ->oowe have 



r 



( J d/A-exp^Qi))^ « exp(- 7 C 2 (r)) 



up to an irrelevant constant, 7 = |^-. Here, C^ir) is the quadratic Casimir operator 
which is, on other hand, the eigenvalue of the lattice Hamiltonian in the strong cou- 
pling approximation fl30l) . Thus we can see the full equivalence between the lattice 
Hamiltonian formulation and lattice in the Euclidean space. 

Two points we would like to point out once more are: 1) the integration over 
space gauge field plays in the Euclidean space a role of the projection onto local gauge 
invariant states; 2) just (and only) this integration or gauge transformation carried out 
at projection ([14]) produces the interactions of the eigenvalues of the Polyakov loops. 

Certainly, all described above is quite known in the context of the lattice gauge 
theories. Reminding this mathematical picture we wanted to make the readers sure 
of two facts: 1) the lattice Hamiltonian formulation is quite reliable approach for 
calculation of the effective action for the eigenvalues of the Polyakov loop; 2) relying 
on apparent gauge invariance of the lattice approximation on the whole we can claim 
that if < A >^ appears it will be a gauge- invariant phenomenon. This question is 
a rather important one and it is, basically, the central point of the loop calculations 
presented in the next chapter. 

Now we are able to consider the problem of constructing of the effective action for 
the eigenvalues of the Polyakov loop. At first we consider pure gluodynamics. The 
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quark contribution is calculated in the chapter 7. The first calculations of the conden- 
sate on the lattice were presented in where the authors carried out the correspond- 
ing Monte-Carlo computations for SU(3) gluodynamics. Although the computations 
were done on the small lattice the unambiguous result was obtained: < A Q >^ in 
the deconfined phase and the condensate falls at the deconfinement phase transition. 
It signifies that Z(N) and SU(N) global symmetries are broken at the same criti- 
cal temperature. Lastly, the similar result comes from |50] where a simulation of the 



pure gauge theory fixing the Landau gauge has been carried out. They have argued 
that Aq develops a nonvanishing expectation value above the deconfinement transition 
temperature and have found the breaking of the colour charge conjugation symmetry. 

Unfortunately, we do not know any other attempts to compute < A > by means 
of the MC-simulations. Here, we are going to explore some known analytical methods 
of statistical mechanics and spin systems. The following material is founded on the 

|] where all omitted (especially technical) details can be found. 



papers |§, 0. II > 

Our starting point is the SU(N) partition function (|4"ID , fl4"2|), with the projection 
operator (| 



We define the effective action in the standard way as [54], |4(| 



S eff (< (3gA >- 



where 



Z 4l = SpK°exp(-PH) 



(55) 
(56) 



K l= Y[dU u {x)U u {x) /n^(^)II n (^)^ n+ (^) 

J X,U J X x,u 

6(N<f>-J2<t>xl (57) 

X 

and N is the number of the lattice sites in the spatial directions. In these terms the 
operator ( fiEf) becomes after substitution (fi x — > + a x 

P = n^°, Px= I rf/i(^)exp(-^(0 + a,) a ). (58) 

x J 

Integration in (|5^) is performed over compact measure with the constraint 

£a x = 0. (59) 

x 

From the last two equations we can deduce that the operator (]58|) projects onto local 
gauge-invariant states with non-zero global colour charge. In a sense it clarifies the 
physical nature of /lo-condensate: it is an imaginary chemical potential for the colour 
charge out of SU(N) group. The simplest case when < A >^ may appear is the 
one of SU(2)/Z(2) gauge group (since elements of center are absent in such a theory 
< A >= in the unbroken phase). Let us consider in brief this theory as a simple 
but non-trivial example which demonstrates apparently how the condensate appears 
in obvious manner. The initial action (|51|) for SU(2)/Z(2) theory includes the adjoint 
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representation for gauge field matrix U n [x). It means that we must use the sum in fl56|) 
only over representations which transform trivially under Z(N) rotations. Performing 
the general steps sketched above and summing over all irreducible representations we 
obtain in the form 

Z su(2)/ Z (2) = r jj d(f)x[s . n jj[e 2(e -7 ) 0x _ 0x+n) _ 02 ( e -7 j 0x + (60) 

T- T n 



X 



where 02 is the Jacobi function. Making use the obvious approximation <p x « <ft + 5a. 
and neglecting all fluctuations 5a x (as the fluctuations are small they can not influence 
the final result because for the gaussian free fields the constraint < 5a x >= is 
automatically fulfilled; of course, it may be not obvious and one should argue this 
point by calculating the corrections; we shall mention this point below for more realistic 
SU (2) group) one finds the effective action to be 

- S eff (<j>) = d\n[l - ^AM ] + (1 - d) In sin 2 0. (61) 

We omitted an irrelevant constant and 7(7) = 1 ex P(~ l^i} + !))• Analyzing this 
expression we conclude that there exists such 7 C that at 7 > j c , 0negO in the point of 
minimum of the effective action (0). This conclusion is quite understandable without 
the approximation we have used here. The Jacobi function in (^) which carries the 
sign "+" does not include the constant part of X . This constant part enters in the 
second Jacobi function with sign " — " . As 7^ lowers the Jacobi function it signifies 
that when 7^ 0, Z$ increases and, thus, leads to deeper minimum of the effective 
action than = 0. Unfortunately, such (almost "strict") proof can be done only in this 
simple example. The more realistic case of SU(2) gauge group has been investigated 
in [p^j . As this case will be presented below and in the more general context we 
restrict ourselves to some phrases here. We used the same approximation as described 
above and calculated the corrections coming from Gaussian integration over 5a x . The 
obtained S e ff was investigated numerically. Our conclusions are similar to the previous 
case. 

Lastly, we are studying QCD (SU (3) gauge group) in details by two independent 
methods of calculations in order to be convinced in the results. First, we suppose that 
all fluctuations are suppressed by coupling constant. Recalling that SU(3) representa- 
tions are labelled by two independent indices / = (Zi, Z 2 ) and eigenvalues of E 2 are the 
quadratic Casimir operator eigenvalues C 2 (l) = |(Z 2 + l\ — hi 2 ) — 1 we can perform a 
summation in the partition function (|56|) over all irreducible representations. In this 
way we can calculate the single-site free energy in the form convenient for the following 
analysis 

- = HF(<f>)) = Hj^] (62) 

where for readability we have introduced the notations 

/,x • 2/01 -02x . 2/ 2 01 + 02^ . 2/202 + 01^ /„ Q> 

fi{4>) = sin ( ) sin ( ) sin ( ) (63) 
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and 



S = exp( 7 ) £[-ie 3 (7/3,0)e 3 (7,2$ Q ) - ^63(7/3,^)63(7,^) + 

^63(7/3, $ Q ) £ e 3 ( 7 /3, ^ - + [e 3 - e 2 ] (64) 

d /3><r 

up to an irrelevant constant. Here 6^(7, 0) are the Jacobi functions and $ Q = (p^—(p a at 
J2 Q ®a = (0 a has been defined after eq.([27|)). Thus < > can be found from the local 
minima of 5" e //G- The S e ffc behaviour has been analyzed numerically at d = 3 and 
in the interval of 7 = [1,2] with high precision. At larger 7 > 7 C (small temperature) 
the minima of ln[F(0)] are located inside of the every triangle plotted on Fig.l at the 
following values: < >= 2nk, < (p 2 >= 2ir(k + |), < </> 3 >= 2ir(k — |) and so on 
(six combinations are possible). Obviously, that this distribution is invariant under the 
Z(3) transformations and we conclude that the system is in the confined phase where 
< TrW >= and < (p >= 0. With 7 decreasing (temperature increasing) at 7 = 1.75 
three secondary minima of the function ln[F(0)] develop and get deeper. At the same 
time the initial minimum of ln[F(0)] gets also deeper but it develops slower and at 7 = 
1.62 it disappears at all. But still earlier at 7^ = 1.73 the secondary minima become 
degenerate with the initial one and the system could undergo the phase transition 
presumably of the first order. In these newly developed minima < <fi >^ which 
signals the global gauge symmetry breakdown and forming temperature dependent 
condensate. It is clear that all these phenomena emerge at the deconfinement phase 
transition since above its critical temperature < TrW >^ in any minimum. 

The second approach to evaluate the condensate for the SU(3) gluodynamics has 
been presented in |^TJ. To understand the picture in more details authors of |51[] have 
considered the first non-trivial term in sum over characters and reparametrized the 
conditions 5{N(j) a — J2 X 0S) m terms of the Polyakov loop. For calculation of the effective 
action Bogolyubov's method of quasiaveraging has been utilized. Numerical studying 
of the effective action up to A 12 (where A is an effective coupling constant which is small 
in the strong coupling region) shows that the condensate appears in the deconfinement 
phase and moreover the picture described above is approximately reproduced in the 
framework of Bogolyubov's method and of redefinition of those quantities over which 
we have to minimize the effective action. In a sense, this method is close to those 
developed in |3(| , [[Hp . The difference is that we used the Polyakov loop to express the 
effective action (though in a special parametrization) from the very beginning. Our 
final result is directly opposite to the one obtained in P0| , ||41|| . It might imply that 
Belyaev's method is not self-consistent (a comprehensive analysis of this method we 
give in the next chapter). 

Generally speaking, the situation could turn out to be more complicated than de- 
scribed above. At high temperature the Polyakov loop develops a non-vanishing expec- 
tation value. In quantum theory, the Polyakov loop, being the function of the random 
variable Aq, is an other random variable. It could mean that to obtain the real dis- 
tribution of contributions to the free energy coming from the Polyakov loop and from 
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Ao-condensate we need to construct a general effective action for both these quantities 
and to find minima of such an effective action. In other case it cannot be excluded 
that either some part of the condensate may be "transfered" into the Polyakov loop or 
vice versa. To verify this idea we carried out the corresponding calculation for SU(2) 
gauge group. Our preliminary results are the following |k|. Implanting 5-function 
5(Ncr — J2 X cos (fix) in ( |57D we considered fundamental SU{2) representation in the sum 
over characters and used the standard mean-field approximation to calculate the effec- 
tive action S e ff(a,<fi). Then, minimizing S e ff(a,(fi) we found that both a and cfi differ 
from zero in the point of minimum of the S e ff(a, cfi) in the deconfinement phase. Hence, 
all the conclusions derived above have been confirmed in this approach. It seems to 
us that in compact theories like lattice gauge theory we should use just this scheme of 
calculations, because invariant integration can, in fact, produce the same < a >^ as 
< (fi >t^ (see (|56|)) but entering with opposite sign and, so, cancel the condensate. It 
may be not the case in the loop expansion method, for at g 2 — > the constraint on a 
in ( |56| ) is automatically fulfilled for the gaussian fluctuations. 

We have, thus, found out that the chromoelectric sector of the lattice gauge theory 
generates Ao-condensate in the deconfinement phase. In all considered approaches the 
obtained picture is essentially the same, perhaps up to irrelevant details. Since the field 
A Q and the relevant Polyakov loop transform non-trivially under the centre subgroup 
transformations the condensate carries out the charges of the centre (as matter fields). 
This leads to quark colour charge screening and vanishing the long-range forces which 
explain, in fact, the deconfinement mechanism. 

Concluding this chapter we would like to make some remarks on Ref . [fj(| ■ The 
authors of the paper have considered a partition function for the eigenvalues of the 
Polyakov loop in the continuum theory. Their consideration is very close to ours de- 
scribed here as they have used the temporal gauge A = with projection onto the 
gauge invariant states. The conclusion of the paper is that there is no real condensation 
at high temperature. The basic assumption conjectured by the authors is the cancel- 
lation of the Vandermonde determinant. We, however, reckon that the real situation 
with the cancellation is a little more complicated than described in The possibil- 
ity of the cancellation of the Vandermonde determinant (group integration measure) 
does exist. It can appear just after integration over spatial gauge field, which one may 
see from ([51]) if we consider a contribution only of spatially longitudinal mode. The 
partition function has a form of the character expansion. The character representation 
can be chosen in such a form that denominator of the sum over characters will include 
group integration measure. And, of course, in this sense the cancellation takes place. 
Let us now return to the Hamiltonian version where we must use a projection operator 
onto gauge invariant states in gauge Aq = 0. It is known for a long time (and we 
stressed this point several times in this chapter) that integration over spatial gauge 
fields is equivalent in a sense to projection onto gauge-invariant states in the gauge 
A = 0. We have shown above that it leads to the same character expansion. In the 
Hamiltonian formulation the character expansion results from projection operator and 
summing over diagonal group operators which commute with the Hamiltonian (the 
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last point was missed in the discussed paper). It means the following: Choosing the 
same representation for SU(N) characters we can formally cancel the determinant. 
But on the other hand, it means that we cancel some part of the projection operator. 
The sum over irreducible SU(N) representations which survives after the cancellation 
can never be presented as the old projection operator (it can be easily demonstrated 
mathematically). The inconsistency of Ref.[|40|[ is the usage of the same projection 
operator after the cancellation. Strictly speaking, this is incorrect. Namely this gives a 
possibility to rewrite the partition function in such a form that the constant part of Aq 
will be only at imaginary unit. Thus, there are two possibilities: either one cancels the 
determinant but after that we must build a new projection operator (if it can exist) or 
we can work with the usual projection operator but the determinant is not canceled in 
this case. In chapter 6 we shall give a confirmation what we sketched above of. The 
same determinant appears not only in the Hamiltonian version but also in the static 
gauge for A Q as the Fadeev-Popov determinant. There is no doubt that cancellation 
takes place and besides not only at one-loop level. But nevertheless the condensate 
falls at high temperature. It may signify only one fact: after the cancellation we have 
no old simple representation for the projection operator. 

Let us consider in the static gauge the part of the Yang-Mills action quadratic in 
spatial fields. Basically, the term at linear power in Aq is nothing but Gauss' law. 
Integration over Aq gives as a result the projection operator onto the states where 
Gauss' law is fulfilled. Let us now make an integration over spatial fields. Since we 
have limited ourselves to quadratic part, this integration can be performed exactly. 
Resulting determinant appears to be the Vandermonde determinant but in (—1/2) 
power as usual. In such a way we come to the obvious cancellation. We verified all 



these results in finite-difference lattice formulation in the same static gauge [56]. Then 



in order to come to the partition function g[J we must look for such an expression for 
the projection operator where Aq enters again at Gauss' law in the action. Only then 
we will be able to separate the constant part of Aq in the exponent. It is impossible, 
which of course can be easily seen from the formulae of this chapter. 

Probably, a more appropriate partition function in the continuum theory just for Aq 
gauge field has been conjectured in [57] where the authors have introduced an invariant 
integration over Aq and postulated an effective action for this variable. This theory is 
able to describe confinement in QCD. We think that theory can exhibit spontaneous 



symmetry breaking by means of Aq field condensation at high temperature [56 



5 < > in loop expansion 



First, a possibility of Ao-condensation had been discussed in the standard loop ex- 
pansion approach to the calculation of the effective action 
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Just in this 

method the most essential questions such as gauge invariance of the condensate, the 
higher loop contributions, the thermodynamics of the Z(iV)-phases have been exam- 
ined. In this chapter we are going to go in a systematic way through the results of the 
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calculations and to discuss them. Our actual calculations will be done for QCD in the 
limit of high temperatures when the coupling constant is small and it is possible to ex- 
pand in this parameter. All calculations are carried out in the background relativistic 
gauge Rf xt which incorporates the results of all relativistic gauges. 

We begin with calculation of the two-loop effective action W^(A , £) of the gluonic 
external field due to the pure gluon contribution [p6fl , [p8] , |27j, |3*I| , |32| (the quark 
contribution will be analyzed in chapter 7). The QCD Lagrangian in the relativistic 
background gauge reads: 

£ = Ik„) 2 + -^(d*q«T + xD»d, x 

+* a ( 7/ A + im)m a + ig^iAl + g^(t c )^ 6 , (65) 

G% = {D B XQ b u-( D u) ah Ql~9f ahc QlQl, 

{D B )ab = S a b ^ + gf a b c A ^ (D ll ) ab = 6 ab d lt + gf abe {Q° tl + A c ri, 
A; = V(^o + ^o), 

where {t c )l are the SU(3) generators, Q a ^ is the quantized field, f abc are the structure 
constants and x, \ are the ghost fields. In what follows it will be convenient to introduce 
the "charged basis" of the gluonic fields: 

^ = ^=(Al± i Al) , 7rJ = 4 , K± = -^=(Al±iAl) , 

Kt = ^{Al±iAl) , V , = Al. (66) 

In this basis the background fields Aq and Aq enter into the momentum space La- 
grangian as constant shifts of the zero momentum components. More details about the 



"charged basis" (p6|) and the Feynman rules are given in |£6 . 

To determine the vacuum value of Aq one should calculate the effective action 
W(A Q ,^) and find its minimum point (A ) min via the minimization procedure. If this 
value occurs to be different from 2|p it means that spontaneous breaking of the gauge 
symmetry happens. 

The effective action W(A ,^) is given as a functional integral over periodic gauge 
and ghost fields and antiperiodic fermion fields (p3|) 

exp[-W(A )VT] = N [ VQV^V^VxD X e' ^ d " x{C+J ^\ (67) 

where £ is the Lagrangian (^), N is a T-independent normalization factor, V is space 
volume and is an external source. The effective action due to gluons up to two-loop 
order in the background gauges has been calculated in Ref. PS|| , f26|j , f27| , plf , p2| , 



34J, | 35| , [33], [55|. Since the result has crucial significance for what follows we have 



calculated it once more. Our result for W(A 0l ^) is the following 

4tt 2 1 3 a 2 3 

W(Ao, OP 4 = + E B ^)) + tEI^V.) + 2B 2 (0)B 2 (a t )} + B 2 ( ai )B 2 (a 2 ) 

6 6U i=i £ i=i 
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+B 2 (a 2 )B 2 (a 3 ) + B 2 {a 3 )B 2 { ai )] + LJ. g p[B 3 (a 1 )[2B 1 (a 1 ) + B 1 (a 2 ) - B^)) 

+ B 3 (a 2 )[2S 1 (o 2 ) + 5 x ( ai ) + S 1 (a 3 )] + B 3 (a 3 ) [2^(03) + 5i(a 2 ) - ^(ai)] (6* 
where the notations have been introduced 



a; = , y - 



7T 7T 



X 1 



01 = -, a 2 = + v%), o 3 = i(-ac + v^y), (69) 



This expression differs from that of ref. |27| in the £- dependent part by the additional 



factor 3/2 . Besides, all signs in the squared brackets are "plus" in p7fl . Our result 
flfi8|) coincides (up to the sign definition) with the result of ref. jJvJ. As is seen from 
flB8|) , the gluon contribution depends on the gauge fixing parameter £. In QCD at high 
temperature the coupling constant is small. So, one can calculate the minimum point 
of W(Ao,£) and the value of the functional in the minimum by an expansion in g 2 . Up 
to the second order we obtain: 

P 4 W gl = -^tt 2 + l -g 2 - J^(3 - 2 awn = J^(3 - 0, = (70) 

Here, the values of x m i n and y m in have been found for the intervals < a^ 2 < 1,-1 < 
a 3 < 0. Five other minima in the (x, y)-plane can be found by means of consequent 
rotations of the coordinate system by the angle |. From ( [70]) it follows that the presence 
of the condensate lowers the action. So, the spontaneous generation of the condensate 
takes place. It is very essential that the depth of all minima is the same. Hence it 
follows that although the condensate breaks both gauge and Z(3) global symmetries, 
the generated phases preserve the rotation symmetry as in the case of Z(A r )-phases. 
The description of the general structure will be done below after evaluating the quark 
contribution. These results are in the full accordance with those from the previous 
chapter. Nevertheless it is to be emphasized that in the loop expansion treatment 
there is no phase boundary whereas in the strong coupling lattice approach there is a 
phase transition to the state with < A >^ 0. 

As is also seen from Q70"D, the vacuum value x m i n as well as the minimum value 
of the functional W(x m i n ) turn out to be dependent on the gauge fixing parameter £. 
This is the ^-dependence problem discussed by many authors |27[| , [ [31H , fl3~2"] , f35|, 



41]. In general, in gauge theories there are two ways of dealing with the problem 
of gauge invariance: l)to use the explicitly gauge invariant formulation from the very 
beginning; 2) to apply the Green function analysis and extract the gauge invariant 
results for the observables via the Nielsen (the Ward type) identities. Besides, some 
hybrid method can be developed as well. As far as the problem of the Ao-condensate is 
concerned, methods of all kind have been applied and different conclusions have been 
extracted. In the previous chapter the explicitly gauge invariant approach to A Q was 
presented. In that scheme < A >^ is obviously a gauge invariant phenomenon. 
Here, we are going to explore the Nilson identity to demonstrate gauge invariance of 
the condensate. 



25 



One of the powerful method of dealing with the ^-dependence problem is to apply 
the Nielsen identities method which first was used in the investigations of spontaneous 
gauge symmetry breaking by radiative corrections f42fl , H58| , |59fl . As is well known, 



in models of the Coleman- Weinberg type both the effective potential V(Q, £) of the 
scalar field $ and its minimum point $ m j n are ^-dependent. A very elegant resolution 
of this problem had been found by Nielsen |42| who proved the gauge invariance of 



this dynamical phenomenon. Namely, he formulated the Ward type identity describing 
an analytical dependence of V A ($,^) on £. So any variation of the potential can be 
compensated by the corresponding variation in (along some characteristic line in 
the ($, £)-plane). More details about these results can be found in [JEJ], ]59fl . In 



the papers [p^l , [pq| this method has been applied to the problem of ^-dependence 
of the gluon condensate. First of all this gives a possibility to check the correctness 
of the calculations. Secondly, due to generality of the approach the relation to other 
methods of calculations can be established. Nielsen's identities of a general form have 
been recently derived by Kobes, Kunstatter and Rebhan For the effective action 
they describe a variation of due to a variation of the gauge fixing term F a ($) 

and in the condensed De Witt notations are given by the expression: 

5W($) = W d 5E j ($). (71) 

Here is a gauge field, superscript " i" includes all discrete and continuous variables, 
denotes a vacuum value of the field, comma after W means the variation derivative with 
respect to corresponding fields and the contraction means integration over continuous 
and summation over discrete variables. The variation 55* describes changing of the 
mean field value due to the special gauge transformations 

S¥ = DiitySQ" (72) 

with the parameter 

5n a = -A^($)5'F /3 ($). (73) 
This parameter is chosen to cancel the change 5S of the classical action 

S g . f .{*) = S{$) + -r^F^)^*) (74) 
due to the variation of the gauge fixing term, 5F a (Q), and the "metric" r] a ^ — > r\ a ^ + 

In Eqs.( |T2"| - [r5| ) -Dq($) are gauge transformation generators, A^($) is the ghost propa- 
gator in the external field S(Q) is the classical action of the gauge fields. In quantum 
theory SE l has to be calculated from the equation |Kj 

5 X l = - < D l a ($)A a J$)5'FP($) > (76) 
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where 

< 0($) >= e- w{ ® J D<$>0(<5>)det[F«(<$>)Dl] exp(-5 g ./. - W d ($)($ - <5) J ) (77) 

and we substituted the current Jj = —Wj. In the background field gauge the gauge 
fixing functions depend on an additional external field $\ So the expression (|76|) should 
be replaced by 

ay = (cj)- 1 ^' (78) 

where the function 

q($) = 5j- < i%(m?m SFP ^;* ) + 

^l^'c*. *)>*=* ( 79 ) 

describes the dependence of the function F a on $. The additional external field must 
be set equal to the vacuum value $ at the last step of the calculations. Besides, the 
latter has to be calculated from the effective action W / ($, $) with $ and $ taken to 
be different. 5j is the Kroneker delta. In QCD one must use the total gauge field 
A" = Ql + b^{b aZ A\ + S a8 A 8 Q ) as the field $\ One must substitute D*, = Df(A + Q) 
as generators and n af3 = j5 ab as the metric tensor. Variations of £ can be realized by 
the variations of the metric 

5r(# = 5'F a = (80) 
In the papers [^lj , [|32]1 the relations (|78"D, ( [79|) have been adopted to the case of the 



background field gauge. To do it the external field $ = A appearing in the gauge fixing 
function F Q ($, $) = (D^(A)Q fJj ) a should be initially assumed to be different from the 
actual background field B a = A a = const. The former must be identified with the 
latter one at the end of the calculations. So the field A in the covariant derivative 
D^(A) can be written as A a ^ = (B + q)^ where are the deviations from the solution 
B®. Taking to be small one may calculate propagators, vertices, etc as series in g° 
and keep only the linear terms. Then the functions b\ n (178|) and Cj (^) are calculated 
by differentiating with respect to q^. On this way the identity ( |7T|) can be written as: 

5W(B) = W, q (B, q) |, =0 [C a JB, q)]~l 6 X ^B, q) \ q=0 (81) 



where 



8x?{B) = -~| < D*{B + Q)Aj(B + q ){D»{B + q)Q;) q=0 >, (82) 
C% = 5« + gf ab %o^ < D a e (B + Q)A e d (B + q)Q; > 9=0 . 



The average values in (ffq), ( p2|) should be calculated in perturbation theory considering 
quantum fluctuations and deviations g" to be small. The internal index "a" in eq. 
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(|ST|) takes now the values a = 3, 8 in accordance with the structure of the background 
field. 

As is seen from eq. (|68|) , the dependence of the effective action on £ appears in the 
order g 2 . So, in the lowest order both the left-hand side (LHS) and the right-hand side 
(RHS) of eq . (jSlD should be of the order g 2 , as well. The covariant derivative D®(A + Q) 
contains the quantum field as a product gQ a ^. So, in the one-loop approximation the 
variation is of the order g. The one-loop effective action W^'(B) has the zero 
order and its derivative with respect to B has the first order in coupling constant. 
Hence it follows that in the lowest order one must use in the identity ( |5I] ) the one- loop 
functions 5x^ and Cg Q to be equal to unity. 

To calculate 5x a (0) it is necessary to take into account the explicit form of gen- 
erators and the fact that only the third and eighth components of external field are to 
be non-zero. Thus, in the one-loop approximation the expression ([82]) reads 



5x (1)a = -\gjf ahc < A M (x - y)Di\B)G%{x - y) > . (83) 
The necessary structure constants in the basis fltjo]) are 



and the background covariant derivatives may be written as follows 

(Z£<#) = D?I% (85) 

where IT* is the column 

n b u = (n + ,n-,n ,k + ,k-,k + ,k-, V )l (86) 

and values of the variables a, b now are: a,b = 7r + , 7r~ , Diagonal elements of D'r 
are the following: 

diagDf = (D«\D- a \d^ D a2 , D~ a2 , D a3 , D~ a3 , M ) (87) 



where Oj stand for the background fields ( |69"D describing the contributions of the isospins 
J, V and U spin subgroups of the SU(3) group, respectively. Using the explicit forms 
of the gluon and ghost field propagators 

G ab ( x - v) = - V f e~ ip{x - y) 

Sab , 



i - 



S db ( x -v) = - V / a p e -W*-v)JL (89) 
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where 



and 



P e ^ = [p° + a e ,p\, p° = — ,n = 0,±l,... (90) 



a 7r± = ±#5 3 , a fc± = ± 9 -{Bl + v^So 8 ), a S± = ±|(-5 3 + v^), a* = dP = (91) 



and substituting ( pq) -(|89]) in eq . (p3[) one obtains: 

^ = if^ [jBl(ai) + ^ l(a2) ~ \ Bl{ - a ^ 

SX 8 = ^y[5i(a 2 ) + SxCas)]^. (92) 

The derivatives of the one-loop parts of W(Ao,£) in eq. (|68D with respect to £> 3 = A 
and B 8 = A equal the expressions: 

dWi _ Agir 



dBl 3/3 3 



[2B 3 ( ai ) + B 3 (a 2 ) - B 3 (a 3 )], 



dW gl 



dBl = ^^[^M + ^ 
By summing up the corresponding products of expressions (92), (93) one obtains the 



gluon contribution to the RHS of eq. fl81|) . The obtained expressions coincide up to the 



sign with the derivative of ( 68f) with respect to £. Thus, the Nielsen identity 



dw 9l _ dwi awj (1) awb (1) _ 

^T" — + W 3 + W 8 " (94) 

is satisfied up to the two- loop order. In eq.(Q) we have denoted C^g = The 
identity (|94|) is just the characteristic equation. So, from the equation it follows that 
along characteristic lines in the (A , £)-plane the effective action is not changing. In 
particular, this is the case for its minimum value (0). In accordance with general 
theory of Nielsen's identity approach this means that the gluon condensation at finite 
temperature is a gauge invariant phenomenon. The fact that identity (|94D holds means 
as well that loop expansion of W(Aq,£) is the self-contained procedure and no other 
^-dependent diagrams should be included in the order g 2 . Only ^-independent terms 
may be added, in principle, to eq . (|68D . From Refs. J42|, |5l|] it also follows that along 



an orbit which passes through the point (£,x m j„ ^ 0) not only the minimum value of 
the effective action W(£, (Ao)mm) but all other observables (particle masses, S'-matrix 
elements, etc) are to be constant as well. This property selects the unique orbit among 
other ones. 

As it was mentioned before, £ dependence of the effective action ( |S8D has called 
in question a possibility of the gluon field condensation because physical phenomena 
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must be gauge independent. So, in a number of papers some gauge invariant methods 
of calculations have been proposed. Historically the first of them has been used by 
Belyaev ([[30]) who introduced a special reparametrization of the background field A 



in terms of the Polyakov loop (18). Applying this procedure to SU{2) gluodynamics 
he came to the conclusion that there is no real condensation at the two-loop level. 
The same result has been derived by this method in SU(N) gluodynamics in Ref.pl 



Thus, it turns out that two methods of calculations give opposite conclusions about 
Ao-condensate (at two-loop level). Let us try to find the origin of the discrepancy. 

First of all let us briefly describe Belyaev's method and consider SU(2) case for 
simplicity (/-spin subgroup of the SU(3) group in (|5"8"|)). The main idea is to define 
the background field x in ( |68D in terms of the Polyakov loop < TrW > with quantum 
fluctuations included. This "classical" or "measured" value, x d , can be calculated via 
< TrW > and in the one-loop approximation the relation between x and x c i has the 
form 

x = x cl + g 2 f(x d ) + o(g 2 ) (95) 

where f{x c i) is a function which has to be found. After the redefinition the effective 
action reads 

W(x) = W(x d ) = W^\x d ) + g^^-W^\x cl )f{x d ) + W®(x d )] + o{g 2 ). (96) 

It has been expected that at least the minimum position in W(x d ) is gauge invariant. 
After calculating the Polyakov loop in the one-loop approximation the following relation 
has been obtained: 

q 2 

x = x d + f^B x (xd/2) £ - f o) (97) 

where £o is an arbitrary fixed number and x = x d when £ = £o- It is the gauge where the 
renormalization of the gluon propagator is absent. Substituting ( p7|) into the effective 
potential and taking into account the explicit form of the Bernoulli polynomials, the 
final result for the two-loop effective action has been found 

^W(x cl ) = vr 2 [-l + ±x d (x cl - 2) 2 ] + <7 2 [1 - ^x 2 cl (x d - 2) 2 ] + o(g 2 ). (98) 

This functional has a minimum at x d = 0, which means the absence of the condensate 
at two-loop level. Now, we analyze this procedure in more details. Actually, it in- 
terpolates between Green's function methods (described above) and completely gauge 
invariant methods like the lattice formulation discussed in the previous chapter. Re- 



ally, in the Refs. |30[ , |41|] as a first step the Green functions are used to calculate the 
effective action. As a second one the reparametrization of the background field in terms 
of < TrW > is performed. The latter may occur to be not consistent. The point is 
that in the loop expansion method, as a rule, the background field A is chosen to be a 
solution of classical field equations and is to be considered as a fixed parameter through 
all calculations. Its vacuum value must be determined via the minimization procedure. 
On the other hand, if one expresses Aq in terms of colourless parameter < W > (but, 
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in principle, other parametrizations are possible) with quantum corrections to be taken 
into account and substitutes it as an argument in the effective action then the variation 
of W^(A ) with respect to < W > may not correspond to the determination of vacuum 
value < A Q > considered as the dominant classical configuration in the initial partition 
function. It seems to us that it would be more natural to express W(Aq) in terms of 
< W > from the very beginning and apply a loop expansion to this functional. This 
scheme was described in the previous chapter and we shown that the condensate does 
appear in such a treatment. 

After these remarks it will be very instructive to describe Belyaev's result in terms 
of variables on characteristics. The basic point here is that the relation ( P7| ) and 
"observable" background fields coincide identically with the characteristic in the 
plane which passes through an arbitrary point (xq = x c i, £o)- I n order to find the value 
of W(x,£) on the characteristic one must substitute eq.(|97|) in eq. (|68D (for the /-spin 
subgroup only in SU(2) case) and then expand in powers of g 2 to order g 2 . After that 
one obtains the equation 

2vr 2 a 2 
(3 4 W charact {x ,£) = —[£ 4 (0) + 2B 4 (x /2)] + y -[B 2 2 {x /2) + 2B 2 {0)B 2 {x Q /2)} 

+ ^g 2 (l - £o)S 3 (V2)Bi(ar /2) (") 

As has been expected, on the orbit the effective action is independent on the parameter 
£ but it is a constant depending on £o,£o- The minimum position and minimum value 

of W charact are 

(Zo(£o))mm = #^(3-£o), (100) 

22 4 

P 4 W charact ((XoUn^ ) = ~ + ~ 7^(3 -£o) 2 - (101) 

15 12 Li)Z7i z 

Taking into account that (xo)min can be identified with x c i one comes to the conclusion 
that (xo)mm in eq.( |101| ) is the gauge invariant "measured" value of the gluon conden- 
sate. This is the exact meaning of the ^-independence of W c haract(xo, £)■ Moreover, in 
this way the idea to express W(Aq,£) in terms of the Polyakov loop can be realized in 
the Nielsen identity approach. As we remarked before, this possibility appears due to 
the fact that characteristics in the (x, £)-plane and the relation of x c i and x are given 



by the same equation (|97|). If one puts £ = 3 in eqs.(p9|), ( |101| ), the effective action 
( j98[) and other results of Refs. |30| , |4l| immediately follow. Hence it is possible to 



conclude that Belyaev's method corresponds to the choice of the orbit which passes 
through the point x c \ = x = 0, £ — 3. This is a special gauge where x c \ 7^ might 
be determined in the three-loop approximation. All other gauges signal < A >^ at 
two- loop level. 

In Ref. ||33|| , the gauge-invariant thermodynamical potential Q(Aq) has been 
calculated step by step in perturbation theory. In this method the minimum position 
of the next order of the potential Cl^ n '(Ao) has to be calculated via the effective action 
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W( n -V(A ,£) of the previous order. In this way the result < Aq >= in two- loop 
approximation has also been obtained. But at the same time the non-trivial minimum 
position for the three-loop thermodynamical potential H^(Aq) was found. Hence, a 
good chance for determination of the gluon condensation at this level has appeared 
in the gauge-invariant method of calculation. Thus, we come to the idea that a final 
determination of Ao-condensate needs a calculation of higher loop contributions. This 
is not a surprise because < A >^ is the two-loop effect and to prove it the first 
quantum corrections should be calculated. This very complicated mathematical task 
has not been solved yet. At the present time in the literature only some partial results 
of the role of the higher loops have been reported. For completeness we shall describe 
them here. 

Generally speaking, when the Green function method is used all the observables 
and, in particular, the higher loop contributions should be calculated with Nielsen's 
identities taken into account. In calculation of W(Aq,£) it may occur that some di- 
agrams of special kind give ^-dependent contributions. Then the only way to check 
the correctness of the approximation scheme is to apply the Nielsen identities. From 
this stand point it is very essential, as we have been convinced before, that the loop 
expansion is a consistent approximation scheme and no other ^-dependent diagrams 
should be added to W^(A ,£) in order g 2 . Only ^-independent contributions can be 
included. Keeping these arguments in mind, let us discuss the results of Refs. [j29f ,|55| 



where the ring diagrams (plasmon diagrams) Wd(A ) have been calculated in the Feyn- 
man gauge (£ = 1). As the main result it has been found that when one considers the 
sum of Wd and the solution < Aq >= in the order g 2 follows. However, from 
the above remarks it is clear that the conclusion contradicts to what Nielsen's identities 
tell us. To understand the situation in details we have calculated the contribution of 
the ring diagrams Wd(A ,^) in the R^ xt background gauge. It was found that this 
contribution is ^-dependent as well as W^(A ,£). Moreover, if we consider the sum 
H/( 2 )(A ,0 + W D (Ao,£) the result < A >= follows again, as in ||29|| . At the same 
time, when one substitutes Wd(A ,^) in the identity it does not fulfilled. So, the 
contribution Wd{Aq,£) is inconsistent. The resolution of this contradiction lies in the 
way of calculation of Wd(A ,^) in [p9fl , |55| . 

As is known [[J, in QCD and QED (with A = 0) the sum of the ring diagrams 
describes the contribution of infrared divergencies to the effective action and results 
in the non-analytic term of order g 3 . In this case the infrared limit is calculated as 
follows: k n=0 — 0, | k \— > ||. The same definition is used in in the case of A ^ 0. 
But actually this may not be the case and an other definition must be considered. In 
|61| the following one has been introduced: k = fc 4 + g/3A — 0, | k |— > 0. Besides, the 



additional contribution resulting from the transversal part of the polarization tensor 
has also been calculated After that the obtained corrections appear to be gauge 



independent and in the limit A — > reproduce the result for the A = case. So, 
just this non-local term should be included to the effective action W(A ,!;). With this 
contribution we have no y4 -condensate elimination. 

Other essential result for understanding the properties of the higher loop effects 
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has been reported in |33| where, in particular, the position of the three-loop thermody- 
namical potential £l^(A ) was determined which is both non-trivial and ^-dependent. 
Hence, the hope to determine < Aq >^ at the three- loop level has obtained a real 
confirmation. Anyway, the calculation of the three-loop contribution to the effective 
action remains of great importance for the final solution of the discussed problem. 



6 Effective Lagrangian approach to Aq condensa- 
tion 

It is well-known that gauge fields are invariant under global Z(iV)-transformations. 
But at finite temperature the Polyakov loop transforms under Z(N) g i subgroup as 
matter fields in the fundamental representation (|l9l). This implies that eigenvalues of 
the Polyakov loop or Aq have the charges of centre of SU(N) like quarks. Besides, the 
Polyakov loop transforms under SU(N) g i group like the matter fields in the adjoint 
representation (P7p. It is a motivation to consider Aq as the Higgs field. Basically, the 
action S G ( 25|) has the Higgs form for (p x defined in eq. fl26|) . These fields are living in the 



lattice sites, they have non-trivial self-interaction and can provide a minimum far away 
from zero. So the action S G is more likely to be the "standard" £77(2) Higgs lattice 
action for the yj-fields. The evidences that at high temperature Aq indeed behaves as 
the Higgs field in the continuum theory and its Lagrangian has the Higgs form have 



been given in|37], |38], |39| . 

In the continuum theory the minimum of the pure gauge action is reached for 
the semiclassical values of the potentials Aq = const, and A°: = 0. If we admit a 
semiclassical value of Aq ^ we immediately find the massive chromomagnetic gluons 
in the effective theory. To obtain the action for the vacuum corresponding to Aq ^ we 
integrate out these massive modes following the philosophy of Appelquist-Carrazone 
decoupling theorem . As a result we obtain an effective Higgs potential for quantum 
fluctuations of field A treated as the Higgs field. This Higgs potential is periodic as a 
function of chromoelectric gauge field with period Due to this periodicity the new 
vacuum of the effective system possesses the periodic symmetry 

ao, A a — A a 4- — 

^0 ^ ^ — ^0 + 

together with the parity symmetry 

Aa ao, _ _ Aa 
A Q ^ A — 

For Aq within the interval of periodicity the effective SU(2) three - dimensional 
theory has been built and it has been shown that the global gauge symmetry is broken 
due to formation of non-zero expectation value of Aq. On the contrary, in the case where 
all the semiclassical potentials vanish the corresponding effective system is completely 
different. We find there massless gluons interacting with Higgs field in a state with 
unbroken symmetry characterized by the effective Higgs potential with global minimum 
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for Aq = displaying only parity symmetry. Transition to the system described earlier 
must be accompanied by a symmetry breaking process. 

Let us now give a sketch of the calculation and show how the global gauge symmetry 
breaking happens. The Fourier expansion with respect to the compactified dimension 
enables us to reduce the four dimensional theory to the effective three dimensional 
field model. The Fourier field components acquire frequency dependent masses and 
therefore it looks reasonable to apply for nonzero frequency modes the Appelquist- 
Carazzone decoupling theorem p2| to calculate an effective theory for static field modes. 
This process generally known as dimensional reduction was described by Appelquist 
and Pisarski|23[] as a possible calculation scheme for high temperature behaviour of 
field theories. Following this approach, an effective static model in three dimensions 
was proposed, reasonably describing the full four dimensional theory in the distance 
range RT = 1 at sufficiently high temperature in the deconfinement phase |38|]. The 
idea of calculation consists in a reduction of a number of degrees of freedom in the 
continual integral definition of the partition function and/or the thermal mean values 
of the observables. The Fourier expansion of the fields with respect to the compacti- 
fied direction was applied which enables us to integrate out the nonstatic modes of the 
fields. The frequency dependent "mass" term prevents the procedure from infrared di- 
vergencies. The ultraviolet divergencies are regularized by zeta function regularization 
method. This procedure results in a nonlocal functional determinant. Expanding the 
determinant we obtain finally the effective action for static modes in three dimensions. 
This expansion is obtained within the background field method for the chromoelectric 
field in the static gauge 

^(x) = < + e(x) (102) 

where £:(x) represents the quantum field fluctuations around the constant field value 
Aq. In the effective theory the chromoelectric field is identified with the Higgs field in 
the adjoint representation^]. The effective action depends explicitly on the static field 
value A through the effective Higgs potential. The mean values of operators calculated 
by this method must be complemented by the information about the value of A in the 
minimum of the effective potential. Depending on the number of degrees of freedom 
chosen we can calculate the different effective systems. To decide which effective system 
for fixed finite temperature corresponds to the full system is the question of comparison 
the minima of the effective potentials corresponding to the different effective actions 
derived from the action of the full system. We suppose that in our approach the 
partition function is not changed by the mathematical operations in the process of the 
calculation of the different effective systems. Therefore the system with lower minimum 
of the effective potential is the most probable effective system corresponding to the full 
theory at fixed temperature. 

Starting from SU(2) gluodynamics we have used the action (|T5| ) expanded up to 
the second order in the chromomagnetic fields. Following the calculation described 
above we come to the effective potential for Aq of the form: 

U eff (X) = -r 4 (^-^[X 4 -2X 3 +X 2 ])-^lnDern gc ° r(X> x2 Xmm ' X) ) (103) 
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/ 2 [(1 - ln(47r) - *(J*f) - (1 -X)} + - 



Z 3 7T 3 




2l 4 



sin 2 (vrX) 



— + 2l 2 | - arctan 

n+X 2 I n + X 



(104) 



Here we have used the new field variable X = gf]A /(2ir)~ 1 , g r is the temperature 
dependent running coupling constant, and 47rl = /3k is the dimensionless momentum. 
The functional determinant in ( |103| ) is calculated by integration over variable \ by 
means of the ^-function regularization scheme. Analyzing the last equation we find 
the following picture. For larger coupling constant (low temperature) the extreme of 
the U e ff(X) is achieved at X = |. This value corresponds to the confinement phase 
where the expectation value of the Polyakov loop is equal to zero. When coupling 
constant is decreasing the minimum value of the potential is removing from the point 
X = | and there will appear two symmetric minima which signal that SU(2) global 
gauge symmetry is broken as in these minima X ^ n where n is any integer. This 
means that in this phase the condensate falls, in full accordance with our previous 
investigations. In the same time we found that the Polyakov loop in this phase differs 
from zero. In this sense this calculations are very similar to those presented in chapter 
4 where the < Aq >^ generation is accompanied by deconfinement phase transition. 
Our approach differs from that of Weiss one |28| by taking into account the quantum 
fluctuations around Aq gauge potential. 

It is interesting to compare the calculations presented above with the similar ones 
done by Oleszczuk in 
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where it has been shown that (Aq) ^ does not appear 
at one-loop level. We find the differences in treating both the chromoelectric part of 
the initial action and massive static modes of the chromomagnetic potentials. In our 
approach the massive static modes of the chromomagnetic potentials were included 
in the effective potential whereas in [39| these modes were missed. As the result the 
Oleszczuk vacuum state contains two states with mass term proportional to (Aq) 2 . It 
implies that the effective potential "wants to correct " the situation and final expression 
for the U e jf does not contain the third order of the background field and possesses the 
minimum for Aq = 0. In this minimum the masses proportional to (Aq) 2 vanish 
and the former massive modes become massless and, therefore, they are a part of the 
vacuum. 

Nevertheless, if we build the semi-classical expansion around the gas of chromo- 
magnetic monopoles||39|| the result (A ) ^ can be recovered. 

As we promised in chapter 4, we are now discussing the problem of the Vander- 
monde determinant cancellation. In the static gauge this determinant appears as the 
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Fadeev-Popov determinant and can be calculated directly without the ghost formal- 
ism. In [[38] it has been shown that the Fadeev-Popov determinant is indeed cancelled 
by contributions coming from the integration over transversal chromomagnetic poten- 
tials. It happens, as was discussed in [JC], even beyond one-loop level. However, it 
does not lead to the conclusion that A$ = in the vacuum [f40[|. In fact the determi- 



nant is cancelled from the effective potential even when Afi ^ fp8|. This fact has 



very simple and beautiful explanation. As we pointed out earlier the integration over 
chromomagnetic gauge field in the present approach is equivalent to projection onto 
gauge-invariant states. It means that if the cancellation takes place when A$ = the 
same must happen also for Aq ^ if the latter is a gauge-invariant phenomenon. From 
this fact we can deduce that the present calculations performed in the static gauge lead 
to the gauge-invariant results since the mentioned cancellation does take place. 



7 Ag-condensate in theory with dynamical quarks 

As it follows from the title, this chapter will be devoted to studying the influence of the 
dynamical quarks on the condensate. We may say that quarks do not change a general 
picture in essential way and both lattice strong coupling approximation with massive 
quarks [[JfJ, []5T| and loop expansion with massless quarks f3"lfl , p^j have shown that 
condensate does not disappear. Nevertheless there are some disagreements between 
these two approaches which should be clarified. Let us begin with calculation of the 
quark contribution on the lattice. 

Two approaches can be developed to solve this task. In |)TJ the hopping parameter 



expansion was utilized to include the effects of massive quarks in Euclidean version 
of the lattice theory. In the quark contribution was obtained in Hamiltonian 



formulation by means of Banks-Ukawa method |25[] and calculation of first not-trivial 
term of high temperature expansion of the fermionic determinant. Since the results 
have been obtained to be essentially the same, we limit ourselves here to elaborate on 
the second approach as it looks as being reliable in more broader range of parameters. 
Lattice Hamiltonian for Kogut-Susskind fermions is of form 

I d _ 

S K _ S = - Vn{x)^{x)U n {x)^{x + n)+m q aY,^{x)^{x), 

x,n=—d % 

V-n = -Vn, U-n{x) = U+(x - n) (105) 

where r] n (x) = ^—\Y 1+X2+ '" +Xn ~ 1 . We remind that we should now use the projection 
operator (^) where 

q a = y( x )\ a ^(x). (106) 

The integration over the quark part of the Hamiltonian yields the factor Z q and for the 
massless quarks (for the sake of simplicity let us take it for a short while) the result 
turns out to read 

Z q = Z° q det(I + 5 ^D? y ) (107) 
q x a a 
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where we have noted 

<£ = l + exp(i^) (108) 

and 



^xy 2 £^,{U%yfiy,%+n U X y5y tX - n ) . (109) 



The diagonal part of the fermionic determinant Z® resulting from the projection oper- 
ator is explicitly given as |25 



2j = n[l + ifefi(&:)]- (HO) 

x 

For non-diagonal part of the determinant the high-temperature expansion can be con- 
structed. The necessary details can be found in p^[ . After calculation we come to 



the following quark contribution restricting ourself to the first non-trivial term in high 
temperature expansion 

- S? ff = InZj + ln( / HdU n (x) exph^Tr^ -U^U + . ...]). (Ill) 

J x,n ° a cr x,n " " 

The dots represent here contributions of the next even terms of logarithmic function 
expansion. All these terms are negative therefore we can regard the maximal values 
of Z® to correspond to the minimal values of the quark effective potential. It is clear 
that maximum of Z® is achieved when < <fi >= 0. This result is very essential as 
it shows the gluon condensate generation is caused by gluonic kinetic energy only 
but not the quark contributions. This point is to be in disagreement with the loop 
expansion method where condensate appears either from gluon or from quark sector 
(see discussion below). 

Further we choose to work with SU (3) gauge theory. As in the pure gluodynamics, 
there are six solutions here, however, due to the explicitly broken Z(N) symmetry they 
are not equivalent. Now the basic maximum is developing at < <p\ >=< 02 >— 0. It 
is clear that evaluating the effect of massive quarks a la Banks- Ukawa i.e. making 
the substitution 

<£^exp(Pm q )+exp(i<p a x ) (112) 

we come to the same conclusion: condensate does not appear at any temperature. 

Performing now the invariant integration in (|111 ) we have up to an irrelevant con- 
stant 

dp 2 [Imn((j))] 2 
16a 2 a [1 + Rett((f))}< 

and then combining fl62|) and ( |113| ) with substitution (|112|) we have analyzed the ef- 
fective potential Sef/G + Sff/ at different values of 7 and m q and found the following 
picture. At larger 7 there is one minimum where condensate is absent. With temper- 
ature increasing (7 decreasing) a phase transition takes place at 7^ which is somewhat 
larger than for pure gluodynamics. Below 7^? the effective action develops two minima 
where < >^ 0, < (p2 >^ 0. It implies that besides spontaneous breaking of the 



S? „ = ln[1 + ^ W ]-^J^L (113 ) 
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global gauge symmetry the condensate generation may lead to spontaneous breaking 
of the colour charge symmetry. This conclusion is in some contradiction to the loop 
expansion method as well. Such a breaking means that in the corresponding minima 
the baryonic number is generated. We elaborate this result in details in the next chap- 
ter. However, at very high temperature when quark mass is vanishing we have found 
that the deepest minimum of the effective action appears to be C-symmetric again 
with < 0i >^ 0, < 02 >= 0. In such a way the disagreement with loop calculations 
can be avoided. Hence, our more essential conclusion is that Ao-condensate falls at 
high temperature in full QCD. < A Q > is non-vanishing because of fluctuations of the 
gluonic kinetic energy and does not present in the quark sector of QCD. 

Now let us present the results of loop calculations in the continuum QCD with 
dynamical quarks. All our notations follow chapter 5 where we discussed the loop 
approach in gluodynamics. The one- loop contribution of massless quarks to the VF(-Ao) 



has been calculated in |L5|. The two- loop quark functional was obtained in Refs. 



35| The Feynman rules and necessary integrals and sums are listed in Appendix A of 



Ref. p4| . Thus, up to two- loop order we have 

3 4 4vr 2 3 
7V / 6 i=i 



X -g 2 { B 2 {a x ) [B 2 ( Cl ) + B 2 {c 2 )\ + B 2 (a 2 ) [B 2 {c x ) + B 2 {c 3 )\ 



+B 2 (a 3 ) [B 2 {c 2 ) + B 2 {c 3 )\ - - £ [ B K c i) - 2B 2 (0)B 2 ( Cl 

^ i=l 

-B 2 {c x )B 2 {c 2 ) - B 2 (c 2 )B 2 (c 3 ) - B 2 (c 3 )B 2 ( Cl ) } 

+ ^y^4 2 { B 1 (a 1 )[B 3 (c 1 )-B 3 (c 2 )} (114) 
+B 1 (a 2 ) [S 3 (ci) - B 3 (c 3 )} + B 1 (a 3 ) [B 3 (c 2 ) - B 3 (c 3 )} } , 

where x, y and are defined in ( |69D and 

!/ 1 !/ 1 

c i = l( x + + 2 ), c 2 = -{-x + -j=y + 2), 

where Nf is a number of quark flavors, Bi(x) is the Bernoulli polynomial of the order 
i defined modulo 1. Calculating the minimum position and minimum value of this 
effective action we obtain 

^W q = -n^N f + 9 ^N f - 9 *N f {3 -^ 
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Here, the minimum values were calculated for intervals 



< di < 1, < a 2 < 1, - 1 < a 3 < 0, 

< ci < 1, < c 2 < 1, < c 3 < 1 (117) 

and, as in the pure gluodynamics, five other minima in (x, y)-plane can be found by 
means of consequent rotations by the angle |. It is interesting to notice that the value of 
the condensate is the same both in the gluon and in the quark contributions. Moreover, 
both fermions and bosons act to lower the action of Aq field. This is rather unusual 
property that the condensate is produced in both sectors of the theory separately. As 
a total result, we have in the two loop approximation 

f3\W Q + W Ql ) min = -7T 2 (— + —Nf) + g 2 (- + —Nf) - g 4 (3 ~? 2 (l + ^) (118) 

HK q-r gljmzn V 45 60 V ^6 72 >' 32vT 2 V <6 ' V ' 

It is worth mentioning that y m in = in the deepest minimum (what means - in the 
vacuum). Hence, it follows that Ao-condensate does not effect a baryon charge at high 
temperature. Now let us elaborate the vacuum structure of gluon condensate which 
results from (|70| ) and ( |118|) . The quarks change the symmetry of the vacuum. As is well 



known, the Bernoulli polynomials are defined modulo 1. Hence, the effective action 
possesses the symmetry in the (x, y) -plane. We display this symmetry in the Fig.l. 
Any dot in the plane can be translated along dashed lines. Vacuum structure consists 
of the hexagonal and triangular elements, which cannot be transformed one to another 
by translations. At the translations the hexagonal elements pass to themselves. So do 
the triangular ones. 

The global minima of the effective potential W q + W g i are marked by dots. Besides, 
in each triangular area there are six local minima which are not depicted. The local 
minima appear only if the quark and gluon contributions are included together. They 
disappear when one considers the gluon contribution only. The local minima depths 
decrease and the global ones increase while the number of quark flavors Nf becomes 
bigger . 

At some Nf the local minima disappear completely. To understand the effect of quark 
mass we also have calculated the contribution ( |1 1 5| ) with m q ^ 0. Unfortunately, this 
calculation is incomplete because of complexity of integrals and the result has been 
obtained for small m q only. It qualitatively looks as follows. The difference between 
the local and global minima decreases with the quark masses increasing. The local 
minima are equal to global ones and we obtain the same result as in the gluodynamics. 
Thus, mainly the light quarks bring up to an appearance of local minima which can 
be identified with metastable phases of the quark-gluon matter. So, one can assume 
that at some intermediate temperature the phase transition from the global minimum 
to the local one takes place. This transition may lead to y m in 7^ in a vacuum. So, it 
can be expected that the chiral phase transition is accompanied by the baryon number 
generation and by spontaneous breakdown of the charge symmetry. 

Now let us check the identity (01) for the quark contribution. Differentiating the 
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one-loop part of W q with respect to B 3 = A and B 8 = A we have 

dW {1) 4vr 

w = ~ww gmci) + Bs(C2) " 2Bs{C3)] ' (119) 

Then by summing up the products of eq. (|92|), (|119|) the RHS of eq. fl8T|) can be 



calculated. The obtained expressions coincide up to the sign with derivative of (|115|) 
with respect to £. Thus the Nielsen identity ( |94] ) is fulfilled for quarks as in the case 
of gluons. Taking into account that non-vanishing v4 -condensate lowers the actions 
W q , W g i and W q + W g i and the above result we see that both basic properties required 
by the Nielsen identities approach are held. In accordance with general theory we come 
to the conclusion that the gluon condensation at finite temperature is a gauge invariant 
phenomenon with quarks included in consideration. 

Thus, both lattice and loop expansion considerations led us to the similar conclusion 
about /Lo-condensate in the full QCD though the effect of the quark contribution to the 
effective action has been found out to be essentially different. This discrepancy does not 
have an explanation at the moment and should be clarified in the future investigations. 
The reason can be the following. In fact, the high temperature expansion used on 
the lattice to calculate a quark contribution is well reliable only at large m q whereas 
the loop calculations have been performed at vanishing quark masses. So it is very 
desirable to carry out both calculations at intermediate values of quark masses. Just 
this point can be the origin of another disagreement concerning a possibility of the 
baryon number generation which will be discussed in the next chapter. 



8 Ag-condensate in hot gauge theories. Some con- 
sequences 

What would be the consequences of such a condensate? We shall present some of them 
and give a brief review of the most important results and related problems. We begin 
our discussion with the problem which was actually one of that difficulties of hot non- 
abelian gauge theories which has impelled to develop the approach to these theories on 
the basis of the possibility of the global gauge symmetry spontaneous breaking. This 
is infrared problem which have been already discussed in different aspects. Now, the 
question is: what may we call a solution of the infrared problem? It has been shown 
in the first papers by Linde that appearing of infrared cut-off even of the order T (the 
condensate appears to be just of this order) cannot save the perturbative expansion 
(in this case all terms of the perturbative expansion starting from the order g G are 
proportional just to this order and therefore all of them give equal in coupling constant 
contributions). From the mathematical point of view the answer is quiet clear - we 
need to construct a method of calculation which would give finite results for expectation 
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values of physical observables in the continuum. One of such methods we discussed 
in chapters 2 and 6. This is the reduction of gauge theories at high temperatures to 
the effective three-dimensional models. It is possible to show that in the course of 
the reduction the static modes of the chromomagnetic potentials acquire mass due to 
formation of the v4 -condensate ||, |39[. Thus, in this method of calculations a hope 



to avoid the infrared divergences by means of the condensate does appear. We think 
that the most essential point in proofing that Ao-condensate can cure the theory of 
the divergences is the proof that just this condensate lead to screening all sources in 
the adjoint representation at high temperatures. The existence of such a screening 
would mean, in our opinion, the solution of the infrared problem. We remind here that 
we should consider spatial correlators to calculate this screening (see our discussion in 
the introduction). Now, we are going to show, omitting all technical details, how this 
screening can appear in the reduction of lattice gauge theories p2| , |14| . 



We present a method much analogous in its idea to the recent development of per- 
turbative expansion resumming conjectured to cure the infrared divergences of 
finite temperature QCD. Actually, the reduction looks like isolation of static contribu- 
tion in the action after Fourier transforming gauge fields A^(x) with further calculation 
of multi-loop corrections over massive non-static modes. This procedure has been al- 
ready performed both in the continuum theory for the Yang-Mills action (see previous 
chapter) and on the lattice for the Wilson action In the latter the reduction 



beyond the perturbative horizon []B5|] , pfifl can be accomplished as the gauge matrix 
U = exp(gaA) expansion is not necessary. 

We develop proper high-temperature expansion and the corresponding reduction 
dealing with the Fourier transformation of compact gauge matrices U^x) rather than 
gauge fields A^x) calculating then a static contribution resulting from such an ex- 
pansion. Two essential points distinguishing this approach from the preceeding ones 
23fl , [j39| , |64fl appear to be the following. Firstly, the static sector generates compact 



dielectric field leading to an effective dielectric theory. Secondly, since in the present 
case Fourier transformation is not simple linear substitution the Jacobian of this trans- 
formation is non-trivial and moreover it generates a mass of dielectric field (which is, in 
fact, the mass of static modes). As to the non-static modes they are massive as before 
2~2" | with the mass proportional to {nT) (in continuum limit). Therefore, in this way 



we could construct quite reliable perturbative expansion for massive modes. All calcu- 
lations of the resulting effective action can be found in , [|l4j] where both corrections 



to the main static contribution and quark contribution have been computed. Here, we 
write down the effective action for static modes restricting to the pure gluodynamics. 

Z= / Y[djj,(a x ) Y[exp(f3' cos<p x cos<p x+n - 2X e cos a x cos a x+n )Z D , (120) 

x x,n 

Zd= I HipfdpMUi)} exp S D (121) 
J i 

f3' = 4(X e /X°) Nt , I = (x,n) is a link. This effective action we call induced dielectric 
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action. In mentioned approximations we have 



where we have laid 



C QStat I QStat I qU /-i 99 \ 



5^ = \nN t Sp'£p(dp)U{dp), (123) 
p 

AT \0 

SL« = —$-E S P mx)*+(x)), (124) 



x.n 



S$* = N t \ e -£Sp(® n (x)V x <S>+(x)V; +n ). (125) 

is zeroth component of gauge field matrices in the static gauge and A m = 

A e = In chosen approximation A ~ 2. The representation $ n (x) = p n {x)U n {x) 

was using in the course of the reduction. 

Comparing formulae (|121]) -( 125 ) with (p!0|)-([T3|) from Introduction one can easily 



deduce that Ql 21 ) is nothing but some kind of LDGT. There exist two differences from 
formulation (13| described in the introduction to be emphasized. First one is that the 
theory ( |121| ) is a compact since the dielectric field obeys the relation < p < 1 unlike 



(|l2l). Another significant difference is the presence of the interaction term S^ at which 
describes interaction between gauge field U n (x), dielectric field p n (x) and Higgs field 
A . 

The mass of the static mode can be calculated from a "naive" effective potential for 
dielectric field. Its classical form can be easily obtained if we neglect the fluctuations 
of the gauge field U n (x) and take the expectation value < A > in S§. The result 
has been put down in eq.(|l4]). It is clear from its form that in the continuum limit 
the mass is proportional to < A > as we need to expand cos^apgAo) when ap tends 
zero (temperature independent constant part should be avoided from the mass after 
adding zero temperature contribution to the free energy). Two facts following from the 
present consideration should be stressed. Firstly, as was conjectured in Introduction, 
two screening mechanisms - dielectric one and that caused by global gauge symmetry 
breaking - can be indeed related as A -condensate supplies the dielectric field with 
gauge-invariant mass. Secondly, we have proved that the mass of the static mode 
appears not only on the level of the standard reduction but beyond the perturbative 
horizon if we take into account all powers of static modes. Now, there is no difficulty 
to prove that spatial adjoint Wilson loop obeys perimeter law at any temperature fl4| . 
Moreover, it is possible to prove that the fundamental Wilson loop displays the area law 



behaviour with string tension to be proportional to the mass of the dielectric field [14 



in full accordance with results of [0, |10[. All of that, according to duality relations, 
mean the screening of the gluon currents and confinement of the "spatial" sources 
being in the fundamental representation. Thus, we can see that A -condensate does 
lead to the screening of the chromomagnetic forces at high temperature. (Of course, 
our proof is not strict and more exact calculation of the adjoint Wilson loop in the 
effective theory with the corrections from non-static modes included is very desirable.) 



42 



The next interesting problem we are going to discuss here is evaluating the heavy 
quark potential in background < A Q >. As is known, heavy quark potential at finite 
temperature is calculated in different phases by using the correlator of the Polyakov 
loops. We find the linearly rising potential in the confinement phase and the Debye 
screened potential in the deconfinement phase where < A >^ 0. 

The heavy quark potential V qq in the finite temperature theory is connected with 
the Polyakov loop correlator as [|17]] 



c 



-^«< x *«> = (L(xi),L(x 2 )> (126) 



In the static gauge we have for the Polyakov loop the simple expression (SU (2) gauge 
group) 

L(x) = cos(i/3s(4) + £(x)), (127) 

i.e. Polyakov loop depends on the chromoelectric field only. This justifies us to use for 
the calculation of the correlator the same method as we have used for the calculation 
of the effective potential in the chapter 6. The effective potential was obtained from 
the relation for the partition function expressed in terms of the effective action 



Z(J3) = J [De]e s °ff (128) 

For the effective action, applying the method described in chapter 6, we have found 
the relation: 



c _ V 



^-^(X 4 -2A- + X 2 ) 
45 3 V ; 



8 J d 3 l s(l)S e £(k,X)e(-\). (129) 



In the last expression e(l) is the Fourier transform of the function /3e(x), S^/(k,X) 
has been put down in ( |105[ ) and all other notations follow chapter 6. 

By the same method one finds for the correlator of two Polyakov loops at the 
distance R the equation: 

(L(0),L(R)) = ie-^ 5 " 1 ^^^ 5 " 1 ^ +cos(27rX)e-# 5 " 1 (°' R )]. (130) 

S -1 means the inverse matrix element with continuum indices given by the relation: 

C-l/ x v)~ f rf3 (ft fc ) 1 r -ik(x-y) /loiN 

The direct application of last formula is not possible on the level of the approxima- 
tions defined in[p8|, because S^J is negative for small k =| k | in the neighborhood of 



X = 1/2. This is not consistent with the definition of physically acceptable action in 
the functional integral. To remove this inconsistency in a natural way we carefully re- 
view the terms abandoned in the previous method of the calculation |38|| with the hope 
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to find some contributions making S^J positive. The calculation is in progress now 
and we would like to stress its significance by pointing out the important consequences. 
Here, we simply modify S^f/ to the positive function S™° d by addition of a small 



positive constant. We rewrite Eq. (131) in the form 



where p = TR. 

The singular structure of the integrand in the complex I plane is defined by the sin- 
gular structure of the S™° d . The function S™° d (l,X) possesses two cuts, (— ioo, —iX), 
(iX,ioo) in the complex I plane. Because the integrands in Eq. (|132|) have no singular- 
ities out of real and imaginary axes, moreover, 



lim I 



SZ d (l,X) ' = 0, (133) 



we change the range of the integration by closing the integration contour. Then we 
can rewrite Eq.( |132| ) after substitution I = ir in the form: 

^(0, R) = J^T tre" r ! m -F ,( :;- X /- (134) 



{27i) 2 pJo \S™ d {lT,X)\ 2 

Applying the mean value theorem, the integral of Eq. (|134|) is estimated as: 

-tiTR 

S-\0,R) = C(T) (135) 

where C(T) is a calculable quantity. The quantity /i is a positive constant from the 
range of integration. For our purposes here it is sufficient to know that p > 0. The 
analytical evaluation of this quantity shows its dependence on the chromoelectric mass. 
Eq.( |135| ) is valid for R ^ 0. The matrix element S ,_1 (0, 0) is infinite and represents 
the quark self energy in the heavy quark potential given in Eq. (|130| ). In what follows 



we shall consider only the interaction energy of the quarks by subtracting S' _1 (0, 0) 
from all relations. We would like to point out the importance of the Eq.( |135| ). All 
qualitative features of the following discussion are based on the particular approximate 
form of the right-hand side of Eq. (|135|) and are independent on details of Eq.( |105| ). 
Investigating the last equations we find the following picture. 

For value of X = | (which corresponds to the vanishing Polyakov loop value (L) = 



0) the heavy quark potential calculated from the Eqs. ( |126|) , ( |130|) , ( p~35|) possesses the 
desired R dependence because for R — » oo the leading term is linear: 

- eract {R) - V q f f ener9y } = -pRT + const - In RT + 0{R- 2 ). (136) 

In the deconfinement phase characterized by (L) ^ the effective potential possesses 
two minima for X ^ 1/2 symmetric regarding to X = 1/2. The heavy quark potential 
in this case is usually calculated by the formula J 



--V qq (R)=ln mo)))2 . (137) 
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Inserting Eqs.( |130|jl35| ) into the last equation we find for R — ► oo the result: 



i n 2 p -fiTR 

-V qq {R) = tg\nX) | C(T) — + 0(iT 2 ). (138) 

In the leading term of the last expression we recognize the Debye screening of the 
colour charge potential in accordance with the lattice results^] as well as analytical 
calculations ||68|| . Thus we can see that appearance of X ^ n (< A >^ 0) leads to the 



Debye screened phase at high temperatures. 

One of the most exciting consequences of v4 -condensate is the spontaneous breaking 
of the colour charge symmetry and, consequently, a generation of the baryonic number. 
If this phenomenon does realize in nature it would be a good possibility both for 
verifying the general theory described in this survey and for searching the new state 
of the strong interacting matter. It should be mentioned at once that the baryonic 
number is generated not at all possible values of the condensate but only when the 
situation described in the end of chapters 4, 5 is realized. In 577(2) gauge theory 
such a generation is not possible at all because the appearance of the condensate 
leads to the symmetric under charge conjugation effective action. In SU(3) theory the 
baryonic number is generated when both < Aq > and < Aq > differ from zero. The 
main statement is: In the presence of a non-vanishing Euclidean (Aq) field the 
generation of the nonzero baryonic number is possible. This statement was 
proved in two cases. 

1) We studied the Hamiltonian formulation of the lattice QCD in strong coupling 
limit []69| where the partition function is of form (the Kogut-Susskind fermions were 
used) ©,(H, (£07]): 

II DKAo) II E exp(- 7 C 2 (0 WA (z)K(A (x + n)) * 

x x,n I 

f J] DU n (x)Det[l + exp(^ Q (x)) + (Ta^ 1 ]T r ]n (x)U n (x)5 x ^ n \ (139) 

x,n n =d 

where 7 = g 2 (2Ta a )~ 1 , Qi is character of the 1-th irreducible representation and the 
high-temperature expansion in the Wilson loops was applied to calculate the fermionic 
determinant. For each colour quantum number the corresponding charge density gen- 
erated by (A ) is given by: 

Pa = p[(A )a,a] = ^(iV)^ ^ ^ ^ ■ (140) 

The chemical potential is introduced into the fermionic determinant. The physical 
barionic density is then obtained as: 

Q = £/>[( A))a,a], (141) 
a 

where a is the colour index. After integrating over gauge fields we get the result: 
Tm W a ( An) 

Pa = ( i + Re ^ ) + (lGT^)- 1 £<A«> + 0(T- 4 ), A a ~ Im W a . (142) 
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From the last equation we can see that baryonic density is proportional to the imaginary 
part of the Polyakov loop. It is easy to conclude that all the higher orders of the (3- 
expansion are also proportional to the imaginary part of the fundamental characters. 
We point out now that (ImW) will be different from zero if (Aq) and (Aq) are non- 
vanishing. 

2). The similar result was obtained in weak coupling regime of the continuum 
QCD for free quarks being in the constant background field Aq |81 |. In that case the 



baryonic number is expressed through different representations of the imaginary part 
of the Polyakov loop: 



l- 



d(A ) a:Ct 

^£ E( . 1) .^ in( „ K) , (143) 



71-2 ^17' " 



Thus we can suppose that in the general well a non-vanishing imaginary part 

of the Polyakov loop yields a non-zero baryonic density even when baryonic chemical 
potential equals zero. It is very important property of (y4 )-condensate which has to 
lead to observable effects in the relativistic heavy ion collisions. 

Now, we would like to draw an attention to the following problem concerning the 
charge broken states. The free energy of SU(3) gauge theory has three global minima 
at the high temperatures. On Fig. 2 we have drawn them schematically for theory with 
dynamical quarks. In the minimum "1" (Tr ImW) = and, therefore, Q — 0. In other 
two minima the charge symmetry is broken, (Tr ImW) ^ and, therefore, Q ^ 0. But 
we do not know in fact whether the C-symmetry broken states are stable or metastable 
and, consequently, short-lived. There are some results concerning this problem. One 
loop calculations ((Aq) = 0) with massless quarks have shown that the states where 
(Tr ImW) = are deeper More precise 2- loop calculations with massless quarks 
(and with (Aq) ^ 0) described earlier lead to the same conclusions. The lattice strong 
coupling evaluations have confirmed the such inference but calculations with massive 



quarks |36| predict the existence of some values of parameters 7 = j3, m(3 (con- 



sidered as independent) preferring the minima "2", "3" to be deeper than minimum 
"1" at 7 — > 0. If it is the case we obtain the following scenario. 

In the chiral symmetric phase at high temperature (m q = 0) the main true mini- 
mum of the free energy will be always real and other two will be only metastable. If in 
the deconfinement phase the quarks have a nonzero mass, the charge broken minima 
become deeper and the baryonic number is nonzero. But at the chiral phase transition 
the system will be back the real minimum again. Basically, in this case the disap- 
pearance of the baryonic number and the transition of the system to the state where 
(Im TrW) = ((Aq) ^ 0, (Aq) = 0) is just the chiral phase transition. To verify this 
picture we can look for the similar phenomena in Z(3) spin system because most of 
features of the Z(3) spin system in d-dimension are similar to the SU(3) ones at finite 
temperature [|T7| . We consider Z(3) theory with external either positive or negative 
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magnetic field, whose partition function is of the form: 



AT-1 n 

z = Y, ex P{ A E M^ qx ■ e^ fe +") +^E C0S (5&)} ( 144 ) 

It is well known that ||72|| : 1) if 77 > 0, the minimum is real, (Im W) = and the phase 
transition of the second order is possible at sufficiently small rj; 

2) if 7] < 0, the minima are complex ones, (Im W) ^ 0, the phase transition of the 
Ising-type takes place to the state with broken symmetry q — > —q. 

What should be concluded from these points? At least at the hopping-parameter 
expansion of fermionic determinant the leading contribution including the Polyakov 
loops has positive sign which definitely leads us to the case rj > 0. Of course, it is 
enough dangerous to do conclusions comparing QCD with Z(3) theory, nevertheless, it 
seems to us that the real minimum will be deeper at any temperature in SU(3) theory 
as well and, hence, the states with nonzero baryonic number are, in fact, metastable. 
In principle, we have not to except other possibility as well. The situation can be much 
more complicated. For example, the discussion in||81|| shows that in QCD slightly above 



the deconfinement phase transition the domains of different "vacua" (with different 
background A ) may exist. The scenario described above cannot be excluded, too, at 
least, until we obtain the connection (or some dependence) between (A ) condensate 
and quark masses. 

Another approach to resolve this problem has been conjectured in [jS2][ : full QCD 
cannot be described by a grand canonical ensemble with respect to triality or quark 
number. The vacuum of QCD with dynamical fermions has triality zero and therefore 
degenerate Z 3 phases and ordered-ordered phase transitions like pure gluonic QCD. 
This means that full QCD does not lead to metastable phases existing up to arbitrary 
high temperatures. 

Certainly, it is very desirable to have more precise both analytical and Monte-Carlo 
calculations concerning the metastable phases with massive and massless quarks to 
clarify this very exciting problem. 

The last question we would like to elaborate here is an inducing of the Chern-Simons 
action in finite temperature gauge theories. We adduce the sufficient conditions of such 
a generation in the background Aq field. 

The Chern-Simons theories are very popular now owing to their unusual and attrac- 
tive properties. The two areas are usually discussed to be relevant to the Chern-Simons 
action: the superconductivity and 4-dimensional field theory at finite temperature. 
There are some well-known examples when the Chern-Simons action is induced by 



fermionic determinant at high temperatures in four-dimensional theory. Redlich|73| 



has studied the model of two left-hand double fermions connected with 577(2) gauge 



field and imaginary chemical potential. In [74] QCD-like theory with axial U(l) charges 



was studied. In both cases the fermionic determinant generates at high temperature 
the Chern-Simons action for the static modes of gauge field. The coefficient in the 
front of the Chern-Simons term is proportional to the introduced chemical potential. 
In connection with it the following question is very interesting: is it possible to induce 
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the Chern-Simons term into finite-temperature QCD and what is the mechanism of it? 
At the first sight this question may seem a bit strange since there are some problems 
which forbid appearance of the Chern-Simons term in the QCD-action: 

1) the Chern-Simons term is defined for 3-dimensional theories; 

2) S c _ s violates the CP-symmetry. 



This part of our survey is founded on the article j75] where all these problems are 
discussed in details. In brief, the S c - S is induced for static gauge modes which are 
described by appropriate 3-dimensional theory. CP-symmetry is broken via certain 
regularization procedure on the lattice. This breakdown vanishes identically in the 
naive continuum limit. The Chern-Simons action is of the form 



S c - S = 8tt 2 k J d 3 x 1(A), 



1(A) = -L £ nmk Sp(A n d m A k + \A n A m A k ). (145) 



Partition function has the standard form 



/ [DA,] exp{- / [dx}L G - T eff (A)} (146) 

J JM 



where T e jf(A) is the result of integration over quark fields. We looked for action ( |145|) 
for smooth potentials of the static gauge modes in T e ff because we do not know how 
it is possible to put down the action ( |145| ) on the lattice in an arbitrary case. We 
used the lattice regularization for the effective action and studied both the Wilson and 
the Kogut-Susskind actions for fermions with different boundary conditions for gauge 
fields. Both the perturbative expansion and the expansion in the Wilson loops of the 
determinant were applied to calculate T e ff(A). The summary of our results concerning 
A -condensate and the Chern-Simons action is as follows: (for detailed calculations see 

El)- 

1. In perturbative expansion a coefficient at the Chern-Simons term is expressed 
through free fermion propagator G, 

K = \l-« ^^mkSpKG-^G^G-^G^G-^G)}. (147) 

Analyzing this expression one can conclude that there is not any perturbative rela- 
tion between A -condensate and the Chern-Simons action (though, exists a possibility 
of k ^ for some values of the Wilson parameter R; see below). Nevertheless, a 
nonperturbative connection does exist. 

2. Either for QCD on the torus due to certain periodic boundary conditions for 
gauge fields 

m = (s 1 )* : a^x + Q = n^x)^ 1 + n^n; 1 ; 

tt^(x + l v )VL v (x) = Vt v (x + l^VL^x) (148) 

or if we take into account nonperturbative configurations out of center of SU(N), 
the fermionic determinant will generate the #-term F^ y F p xe^ vp \. We suppose that 
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the vortex configurations of gauge field, existing on the lattice, can provide the same 
minimum of the effective action as fields = [76] : 



Sp Z{dp)V(dp) = 1, Z{dp) G Z(N), 
V(dp) = V n (x)V m (x + n)V+(x + n)V+(x), 
V m (x) e SU(N)/Z(N). 



(149) 



The #-term is reduced to the Chern-Simons action at high temperatures and when 
A condensate falls. The coefficient at the action is k ~ anf3(A ) for the theory on the 
torus and n ~ \/3 (A ) F nm (@) for the second case. Here, n is the winding number of 
matrices Q and @ is a vortex potential. 

3. In theory with the SU(3) gauge group and nonperturbative dielectric vacuum 
S c _ s is generated when expectation value of the imaginary part of the Polyakov loop 
differs from zero. Nontrivial dielectric vacuum appears in an effective three-dimensional 
theory with the group having nontrivial center f62"fl , ||i~4"f (see also our discussion above). 
In the last both cases the configurations from the center of SU(N) are necessary for 
generation of 5 C _ S . So it is very plausible that such a mechanism is impossible in the 
standard continuum QCD where the variables from the center of SU(N) are absent. 
Besides, there is another necessary condition of the appearing of 5* c _ s here. In any of 
these cases QCD has to belong to definite class of universality. It means the following. 
Choosing the lattice fermionic action we had to solve the problem of doubling the 
fermion degrees of freedom. We should use either the Wilson action [I77J 



Sw = -a 3 ^[^(x)f/ M (x)GR-7 M )^(x + cv) +V(x)U+(x - a ll )(R + 'y lt )V(x - a,*)] 

x) (150) 



+ a i m^2^(x)^(x)-da 3 J2^(x)R^\ 



or the Kogut-Susskind action fl4"5 
Sk-s 



7, Yl Vn(x)ty(x)U n (x)fy(x + n) + ma^ ^(x)^(x), 

x,n=-d x 



rj-n = -Vn, U- n (x) = U+(x- n) 
which allow us to avoid undesirable degrees of freedom. The most general form of the 
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Wilson parameter R is |L2, 

R = s exp{%6^ 5 )T, -< s < 1, < Q < vr, T 2 = 1, 



where T is some matrix from the fermionic space. One of the //-symbols is |79 
Vn (x) = (_i)*i+*2+~+*«-i exp{nr[k{x + n) - k(x)}) 



(152) 



(153) 



where k is any integer. At the level of free quantum theory we cannot choose a unique 
form of these quantities. More than that, it is just valid at any choice of R and r\: 



49 



- naive continuum limit is(6 l , s, T, /c)-independent and chiral symmetry is restored; 

- fermionic propagator and lattice Feynman's rules coincide in the continuum limit 
with standard Feynman's rules; 

- the property of the positivity is fulfilled. 

Nevertheless, the theories with different R can belong to different classes of uni- 
versality since they can define different quantum continuum theories. Examples of the 
such nonuniversality can be found in [jSOfl . In the only case if 77-symbols or parameter 
R satisfy some conditions, for example, 

Vn{x)i] m (x + n)r] k (x + n + m)rji(x + n + m + k) = e nm ki (154) 

the Chern-Simons action will be generated in QCD. Fortunately this property does 
not depend on lattice regularization P0 |. On the other hand it is very difficult task 



to prove that the S c - S is unique nonuniversal contribution to effective action in finite- 
temperature QCD though it seems very plausible (see, for instance, The parame- 
ters R and 77 as in (|154|) violate the CP-symmetry. As we could see above, this violation 



disappears in the continuum theory both in the naive limit and in the free quantum 
theory. The appearance of S c _ s is the only remnant of this violation which survives in 
the continuum quantum limit. There is only one theoretical reason for the choice of 
R or 7], the property of the positivity ||78|| . But this property is satisfied for given R 
and r\. So we may formulate the question: is it possible that QCD belongs to the class 
of universality in which the Chern-Simons action exists? The answer is left unclear 
so far. Certainly, the presence of S c - S has to lead to interesting phenomena at high 
temperatures because S c - S can change a statistics of a matter fields connected with it. 
Thus, appearance of the Chern-Simons term in the QCD action at high temperature 
is approved by the following complex of the circumstances: 

a) appearance of A -condensate; 

b) existing of nonperturbative vacuum which is formed by vortex potentials from 
the center of the gauge group; 

c) problem of universality. 

To finish, we would like to notice that appearing of S c - S could lead to solving of 
infrared problem as well, since the Chern-Simons action generates magnetic mass just 
for the infrared dangerous static modes. 



9 Discussion and Summary. Unsolved problems 

In the present survey we have considered the mechanism of the spontaneous breaking 
of the global gauge symmetry caused by the condensation of the gluon field at high 
temperatures. We have discussed in details the most important approximations for 
calculation of Aq condensate both in the continuum theory and within lattice gauge 
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models. We understand that it is almost impossible to give a full review of all ques- 
tions related to A condensation because breaking SU(N) symmetry is certainly strong 
phenomenon which has to reflect itself on many aspects of high temperature QCD. 

As an essence of the above analysis we would like to stress once more that in all 
considered approaches the gluon field condensation at T 7^ has been determined and 
proved to be a gauge invariant phenomenon. The nice and important fact that in all 
approaches we have discussed here the gauge invariance of the condensate is manifested 
by the different methods makes us sure that the condensation may indeed be realized 
in the nature. At the same time, a number of discrepancies in results obtained by 
different methods (mainly as the role of quarks is concerned) have been found. These 
points need to be investigated separately in order to have more reliable results. 

It is natural that in the every method of calculation there are own problems which 
should be solved in the future. Generally speaking, they can be divided in three 
categories dependently on their importance and theoretical and practical significance. 
The main problem now is to derive a conclusion on A Q condensation with higher loop 
contributions to be included. This is a complicated mathematical task which should 
be investigated on the base of Nielsen's identities in order to obtain a correct gauge 
invariant result. Then, if < Aq >^ will be determined finally, the problem of 
constructing the consistent theory with < A >^ should be solved. This point 
should be investigated in details because just here a number of theoretical problems 
remain unsolved yet. Most essential of them is the construction of partition function 
and calculation of various observables with < Aq >^ to be taken into account. As 
we have seen, there are six minima with the same depth in the (x, y)-plane and Z(3) 
symmetry is broken. This situation, but in the case when Z(N) is preserved, has been 
discussed in the paper |81 |. It was shown that an extension of standard description 



of thermal equilibrium is needed when discrete degenerate phases are to be included. 
However, how does it work in the case of spontaneous breaking of SU(N) has not been 
investigated yet. In this line the investigations of infrared problem of gauge fields at 
T ^ 0, the gluon magnetic mass, asymptotic behaviour of running coupling constant 
are of great importance, too. Having these questions answered one obtains a consistent 
theory with < A >^ 0. After that a phenomenology of finite temperature QCD will 
be more transparent. And only this third of mentioned points may give a possibility 
of a deriving of a final conclusion on gauge field condensate. Monte-Carlo simulations 
and analytical evaluations on the lattice making use some correlation inequalities could 
help in this problem as well. 

At last, the very interesting and exciting problem is to find a dependence of A - 
condensate on quark mass in order to clarify the possible connection between < A > 
and the order parameter of the chiral phase transition. One may hope that such a 
connection really exists since the Ao-condensate can drastically change the structure of 
the quark propagator. On this way we may also hope to find a solution of the problem 
of the baryonic number generation. 

O.A.B. thanks Institute of Physics of Slovak Academy of Sciences for their kind 
hospitality where the main part of this review has been prepared. The authors are 
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